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ABSTRACT 


The  performance  of  two  novel  interferometer  antenna  architectures  for  high- 
resolution,  wideband  direction  finding  are  investigated.  The  first  configuration  incorporates 
a  Symmetrical  Number  System  (SNS)  encoding  of  the  interferometer  amplitude  response 
(symmetrical  folding  waveform).  The  second  configuration  incorporates  a  Residue  Number 
System  (RNS)  encoding  of  the  interferometer  phase  response  (saw-tooth  waveform).  The 
residue  architectures  serve  as  a  source  for  resolution  enhancement  in  an  interferometer  array 
by  decomposing  the  analog  spatial  filtering  operation  into  a  number  of  parallel  sub¬ 
operations  (moduli)  that  are  of  smaller  computational  complexity.  Each  sub-operation  only 
requires  a  precision  in  accordance  with  the  size  of  the  modulus.  A  much  higher  resolution 
is  achieved  after  the  N  moduli  are  used  and  the  results  of  these  low  precision  sub¬ 
operations  are  recombined.  A  four-element,  3  channel  array  using  the  moduli  set  mj  =  3 , 
m2  =  4  and  m3  =  5  was  constructed  in  a  ground  plane  using  rectangular  waveguide 
elements  with  a  center  frequency  of  8.5  GHz.  Experimental  results  are  compared  with  the 
simulation  results  to  demonstrate  the  advantages  of  this  approach.  The  frequency  response 
of  the  RNS  array  is  investigated  numerically.  To  correct  the  quantization  errors  due  to  any 
frequency  offset,  a  fast  correction  algorithm  is  derived  and  is  shown  to  have  excellent 
results  over  a  wide  bandwidth. 
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I.  INTRODUCTION 


A.  OVERVIEW  OF  DIRECTION  FINDING  (DF)  SYSTEMS 

Direction  Finding  (DF)  systems  provide  an  emitter’s  bearing,  or  angle  of  arrival  that 
can  be  used  as  an  invariant  sorting  parameter  in  the  deinterleaving  of  agile  pulsed  radar 
signals  and  in  separating  closely  spaced  communications  emitters.  In  addition,  the 
conservation  of  jamming  power  in  power-managed  Electronic  Attack  (EA)  systems 
depends  on  the  ability  of  the  associated  Electronic  Warfare  Support  (ES)  system  to  measure 
the  direction  to  the  victim  emitter.  Direction-of-arrival  measurement  systems  fall  into  three 
main  categories:  scanning  beam,  simultaneous-multiple-beam,  and  interferometer 
techniques.  The  mechanically  scanning  beam  requires  only  a  single  antenna.  However,  this 
technique  has  the  disadvantage  of  a  low  probability  of  intercept.  The  simultaneous- 
multiple-beam  system  uses  an  antenna  or  several  antennas,  forming  a  number  of 
simultaneous  beams,  thereby  retaining  the  high  sensitivity  of  the  scanning  beam  approach 
while  providing  fast  response.  However,  it  requires  many  parallel  receiving  channels,  each 
with  full  frequency  coverage.  Interferometer  systems,  based  on  phase-comparison 
techniques,  have  the  advantage  of  a  fast  response  but  generally  use  wide  coverage  antennas 
which  result  in  low  sensitivity.  In  addition,  they  require  relatively  complex  microwave 
circuitry  that  must  maintain  a  precise  phase  match  over  a  wide  frequency  band  under 
extreme  environmental  conditions.  When  high  accuracy  is  required  (on  the  order  of  0. 1  to  1 
degree),  wide  baseline  interferometers  are  used  with  ambiguity  resolving  circuitry.  [Ref.  1] 

The  two  primary  interferometer  techniques  used  for  DF  are  the  amplitude- 
comparison  method  and  the  phase-comparison  method.  Generally  the  phase-comparison 
method  has  the  advantage  of  greater  accuracy,  but  the  amplitude-comparison  method  is 
widely  used  due  to  its  lower  complexity  and  cost.  The  basic  interferometer  geometry  is 
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illustrated  in  Figure  1.1  whereby  a  plane  wave  arriving  at  an  angle  0  is  received  by  one 
antenna  element  earlier  than  the  other  due  to  the  difference  in  path  length.  The  time 
difference  can  be  expressed  as  a  phase  difference 

A(t)  =  — r— -sin©  (1-1) 

A 

where  d  is  the  spacing  between  antennas,  and  A  is  the  wavelength.  For  a  two-element 
interferometer,  the  spacing  d  must  be  A/2  or  less  to  provide  unambiguous,  or  single  lobe 
+90°  coverage. 


Figure  1.1:  Phase  Interferometer  Principle. 
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The  spacing  between  antennas  must  be  established  for  the  highest  frequency  to  be 
received.  Interferometers  employing  multiple  antenna  elements  are  called  multiple  baseline 
interferometers.  Usually  the  receiver  consists  of  a  reference  antenna  and  a  series  of 
companion  antennas.  The  spacing  between  the  reference  element  and  the  first  companion 
antenna  is  A/2;  other  secondary  elements  are  placed  to  form  pairs  separated  by  1,  2,  4,  and 
8  wavelengths.  This  multiple  baseline  interferometer  provides  a  binary  angle-of-arrival 
measurement  where  each  bit  of  the  measurement  supplies  a  more  accurate  estimate  of  the 
emitter’s  angle  of  arrival.  [Ref.  1] 

Another  kind  of  multiple  baseline  interferometer  is  the  nonharmonic.  In 
nonharmonic  interferometers,  no  pair  of  elements  provides  completely  unambiguous  angle- 
of-amval  information  over  the  complete  field  of  view.  For  example,  the  initial  spacing  in 
the  nonharmonic  interferometer  might  be  A,  while  the  next  companion  element  spacing  is  3 
A/2.  Ambiguities  are  resolved  by  truth  tables,  and  hence  the  accuracy  is  set  by  the  spacing 
of  the  widest  baseline  antenna  pair.  Nonharmonic  interferometers  have  been  implemented 
over  nine-to-one  bandwidths  (2  to  18  GHz)  with  rms  accuracy  from  0.1  to  1  degree  and 
without  ambiguities  over  ±90°.  The  principal  advantage  of  the  nonharmonic  interferometer 
over  the  harmonic  interferometer  is  the  increased  bandwidth  for  unambiguous  coverage 
[Ref.  2]. 

B.  PRINCIPAL  CONTRIBUTIONS 

This  thesis  investigates  the  performance  of  two  new  interferometer  antenna 
architectures  for  high-resolution  wideband  direction  finding.  The  first  configuration 
incorporates  a  symmetrical  number  system  (SNS)  encoding  of  the  interferometer  amplitude 
response  (symmetrical  folding  waveform).  The  SNS  can  serve  as  a  source  for  resolution 
enhancement  in  an  interferometer  by  decomposing  the  analog  spatial  filtering  operation  into 
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a  number  of  parallel  sub-operations  (moduli)  that  are  of  smaller  computational  complexity. 
Each  sub-operation  symmetrically  folds  the  antenna’s  amplitude  response  with  folding 
period  equal  to  twice  the  modulus.  Thus,  each  sub-operation  only  requires  a  precision  in 
accordance  with  that  modulus.  A  much  higher  resolution  is  achieved  after  the  N  different 
SNS  moduli  are  used  and  the  results  of  these  low  precision  sub-operations  are  recombined. 
Simulation  results  are  used  to  numerically  evaluate  the  performance  of  a  4-element  (3 
channel)  array  using  moduli  0^=3,  1x^=4  and  m3=5.  The  antenna  was  constructed  in  a 
ground  plane  using  rectangular  waveguide  elements  with  a  center  frequency  of  8.5  GHz. 
Experimental  (anechoic  chamber)  results  were  obtained  and  compared  with  the  simulation 
results  to  demonstrate  the  feasibility  of  the  approach.  The  second  antenna  configuration 
incorporates  a  residue  number  system  (RNS)  encoding  of  the  interferometer  phase 
response  (sawtooth  folding  waveform).  The  phase  response  folds  with  folding  period 
equal  to  the  modulus.  Simulation  results  indicate  that  die  RNS  architecture  antenna 
provides  a  more  robust  solution  to  high-resolution  DF.  The  frequency  response  of  the  RNS 
array  is  investigated  numerically.  Quantization  errors  in  the  resolved  direction  of  arrival 
when  the  intercepted  frequency  is  other  than  the  center  frequency  are  demonstrated.  A  fast 
correction  algorithm  that  uses  the  known  intercepted  frequency  (e.g.,  from  an  intercept 
receiver)  is  derived.  The  performance  of  the  algorithm  is  quantified  and  shown  to  have 
excellent  results  over  a  large  bandwidth. 

C.  THESIS  OUTLINE 

Chapter  II  details  the  mathematical  development  of  the  linear  antenna  array.  A  two 
element  interferometer  is  shown  as  a  special  case.  The  array  factor  associated  with  the 
linear  array  is  derived.  Chapter  HI,  is  devoted  to  the  encoding  of  amplitude  response  with 
the  SNS.  First  the  SNS  encoding  procedure  is  briefly  introduced.  The  SNS  antenna 
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architecture  is  then  derived  with  simulation  results  shown  to  demonstrate  the  array 
performance.  Construction  of  the  antenna  hardware  is  described  and  experimental  results 
are  shown.  Chapter  IV  examines  encoding  the  phase  response  with  the  RNS.  The  RNS 
encoding  procedure  is  briefly  described  and  simulation  results  are  shown.  In  Chapter  V  the 
antenna  frequency  response,  antenna  size  and  respective  accuracy  are  investigated  for  the 
RNS  antenna  design.  The  correction  algorithm  is  derived  using  the  intercepted  frequency. 
Chapter  VI  presents  a  summary,  some  concluding  remarks  and  recommendations.  The 
Appendix  provides  the  Matlab™  code  that  simulates  the  output  of  both  the  SNS  and  RNS 
antenna  for  direction  finding,  as  well  as  the  calculations  of  the  quantization  errors. 
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II.  LINEAR  ARRAY  ANTENNA  THEORY 


A.  OVERVIEW 

Several  radiation  elements  can  be  arranged  in  space  and  interconnected  to  produce  a 
directional  radiation  pattern.  Such  a  configuration  of  multiple  radiating  elements  is  referred 
to  as  an  array  antenna,  or  simply  an  array.  In  an  array  of  identical  elements,  the  radiation 
pattern  is  determined  by  the  type  of  individual  elements  used,  their  orientations,  their 
positions  in  space,  and  the  amplitude  and  phase  of  the  current  feeding  them.  To  provide 
specialized  radiation  patterns,  it  is  necessary  that  the  fields  from  the  elements  of  the  array 
interfere  constructively  (add)  in  the  desired  signal  directions  and  interfere  destructively 
(cancel  each  other)  in  the  remaining  space.  Ideally  this  can  be  accomplished,  but  practically 
it  is  only  approached.  Arrays  can  be  found  in  many  geometrical  configurations.  The  most 
elementary  is  that  of  a  linear  array  in  which  the  array  element  centers  lie  along  a  straight 
tine.  The  elements  may  be  equally  or  unequally  spaced.  When  the  array  element  centers  are 
located  in  a  plane,  it  is  said  to  be  a  planar  array.  Examples  of  planar  arrays  are  circular  and 
rectangular  arrays  in  which  the  element  centers  are  disposed  on  a  circle,  or  are  contained 
within  a  rectangular  area,  respectively.  Another  class  of  arrays  is  that  of  conformal  arrays. 
In  this  application  the  array  element  locations  must  conform  to  some  nonplanar  surface 
such  as  that  found  on  an  aircraft  or  missile.  [Ref.  3] 

The  simplest  and  one  of  the  most  practical  arrays  is  formed  by  uniformly 
distributing  the  elements  along  a  line.  This  arrangement  is  the  conventional  linear  array 
configuration.  Considering  that  all  of  the  array  elements  are  identical,  the  principle  of 
pattern  multiplication  applies.  This  is  stated  as:  “The  electrical  field  pattern  of  an  array 
consisting  of  similar  elements  is  given  by  the  product  of  the  pattern  of  one  of  the  elements 
(the  element  factor  (EF))  and  the  pattern  of  an  array  of  isotropic  point  sources  with  the 
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same  locations,  relative  amplitudes  and  phases  as  the  original  array  (the  array  factor  (AF))” 
[Ref.  4].  That  is, 


E(e,<j>)  =  EF(e,<t>)-AF(0)  (2-1) 

where  E(0,  <|))  is  the  total  radiated  (or  received)  electric  field  intensity  in  the  direction  (0 , 

<|>). 

The  array  factor  may  be  thought  of  as  a  sum  of  isotropic  point  sources  located  at  the 
center  of  each  array  element.  The  array  factor  is  that  factor  of  the  radiation  pattern  which  is 
found  from  the  element  currents  (amplitudes  and  phases)  and  their  locations.  On  the  other 
hand,  the  element  pattern  is  that  factor  of  the  radiation  pattern  determined  by  the  individual 
properties  of  an  element  (its  current  distribution  and  orientation  in  space). 

B.  LINEAR  ARRAY  ANTENNA  BEAMFORMING 

A  typical  linear  array  composed  of  Ne  similar  element  antennas  is  shown  in  Figure 
2.1.  The  output  of  each  array  element  may  be  controlled  in  amplitude  and  phase  as 
indicated  by  the  phase  shifters  and  attenuators.  In  addition  to  the  amplitude  and  phase 
control  imposed  on  each  element  represented  by  In  (n  =  0,  1,  2,...,  Ne-1),  there  is  relative 
phase  shift  between  the  waves  arriving  at  the  antenna  elements  due  to  their  positions  in 
space  and  the  angle  of  arrival  of  the  wave.  Furthermore,  the  pattern  of  each  element  leads  to 
a  response  that  varies  with  the  angle  of  arrival  of  the  incoming  plane  wave.  If  it  is  the  same 
for  each  element,  the  principle  of  pattern  multiplication  allows  us  to  consider  it  as  one 
multiplicative  factor  in  the  total  pattern.  The  remaining  angular  dependence  of  the  pattern  is 
called  the  array  factor,  and  it  is  determined  only  by  the  element  positions  and  their 
amplitudes  and  phases  represented  by  In.  The  array  in  Figure  2. 1  is  shown  as  a  receiving 
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array;  however,  the  receiving  pattern  is  the  same  as  the  transmitting  pattern  by  reciprocity 
[Ref.  4]. 


Figure  2.1:  A  Typical  Linear  Array.  The  Output  Currents  are  Summed  in  a  Summing 

Device  Before  Entering  the  Receiver. 

Since  the  pattern  of  an  array  obtained  by  neglecting  the  patterns  of  the  individual 
elements  of  the  array  is  the  array  factor,  then  if  we  replace  each  element  of  the  array  with  an 
isotropic  point  source,  the  resulting  pattern  is  the  array  factor.  The  isotropic  point  source  is 
an  idealized  antenna  occupying  a  point  in  space  and  radiating  uniformly  in  all  directions. 
The  radiation  fields  of  a  point  source  at  the  origin  of  a  spherical  coordinate  system  are 
proportional  to  I0e  *747tr  where  I0  is  the  current  of  the  point  source,  k  is  the  wave-number 
equal  to  2k/X,  and  r  the  radius  of  the  sphere  [Ref.  4].  The  far-field  pattern  is  obtained  from 
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the  angular  dependence  of  the  radiation  fields.  Thus  the  pattern  of  a  point  source  is  constant 
and  given  by 


AF  =  I0.  (2-2) 

As  elements  are  added  to  the  array,  each  one  with  a  different  current,  it  is  necessary  to 
account  for  their  relative  field  strengths  as  determined  by  their  element  currents. 

The  array  factor  for  the  array  in  Figure  2.1  is  derived  using  the  array  in  Figure  2.2 
which  has  point  sources  for  the  array  elements  in  place  of  the  actual  elements. 


Figure  2.2:  Equivalent  Configuration  of  the  Array  in  Figure  2.1  for  Determining  the  Array 
Factor.  The  Elements  of  the  Array  are  Replaced  by  Isotropic  Point  Sources. 
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The  array  factor  for  this  receiving  array  is  then  the  sum  of  the  point  source  responses 
{eJXK\ eJVl , e)VNe~' }  each  weighted  by  the  amplitude  and  phase  shift  {I0,  Ij,...,IN  _J 
introduced  in  the  transmission  line  connected  to  each  element.  Here  \|/0,\|/1,...,\|/Ne_1  are  the 
phases  of  an  incoming  plane  wave  at  the  elements  locations  designated  0,l,...,Ne  —  1. 
Usually  these  phases  are  relative  to  the  coordinate  origin;  that  is,  the  phase  of  the  wave 
arriving  at  nth  element  leads  the  phase  of  the  wave  arriving  at  the  origin  by  \|tN  _v  [Ref.  4]. 
Thus,  the  array  factor  of  the  array  shown  in  Figure  2.2  is  given  by 

AF  =  I0ej't'°+I1ejVl +-  +  INe_iejVN-1.  (2-3) 

Now,  consider  a  linear  array  of  Ne  elements  as  shown  in  Figure  2.3. 


Figure  2.3:  Equally  Spaced  Linear  Array  of  Isotropic  Point  Sources. 


The  elements  are  equally  spaced  with  a  distance  between  elements  equal  to  d.  The  angle  0  is 
that  of  an  incoming  plane  wave  relative  to  the  axis  of  the  array.  The  isotropic  sources 
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respond  equally  in  all  directions,  but  when  their  outputs  are  added  together  (each  weighted 
according  to  In)  a  directional  response  is  obtained.  The  phase  of  the  wave  arriving  at  the 
origin  (n=0)  is  arbitrarily  set  to  zero,  y0=0.  The  incoming  wave  at  element  1  arrives  before 
those  at  the  origin  since  the  distance  is  shorter  by  an  amount  dsin0.  The  corresponding 
phase  leads  at  element  1  relative  to  those  at  element  0  is  kdsinG.  This  process  continues  and 
Equation  (2-3)  becomes 

N.-l 

AF  =  I0  +Iiejkdsine  +I2ejk2dsine  +  •••  =  2X  ejkndsine.  (2-4) 

n=0 

Now  consider  the  array  to  be  transmitting.  If  the  current  has  a  linear  phase 
progression  (i.e.,  the  relative  phase  between  adjacent  elements  is  the  same),  we  separate 
out  this  phase  explicitly  as 


In=Anejna 


(2-5) 


where  the  n+1  element  leads  the  nth  element  in  phase  by  a.  Then  Equation  (2-4)  becomes 


N,-l 


AF=lA,e 


jn(kdsin0+a) 


n=0 


Define  y  as 


y  =  kd  sin0  +  a 


(2-6) 


(2-7) 


then 


N.-l 

AF=  XA.e*”'. 

n=0 


(2-8) 
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A  very  important  special  case  of  equally  spaced  linear  arrays  is  that  of  the  uniformly  excited 
array.  This  is  an  array  whose  element  current  amplitudes  are  identical,  so  that 

Ao  —  Ai  —  A2  =  •••  =  An^_j.  (2-9) 

The  array  factor  from  Equation  (2-8)  is  then 

N.-1 

AF  =  A0£e^=A0[  l  +  ejv  +  ej2v  H - 1- ej(Nc-1)'t'] .  (2-10) 


A  few  short  steps  are  required  to  sum  this  geometric  series.  First  multiply  the  last  equation 
by  eJV  to  obtain 


AFejv  =  A0  [ejv  +  ei2'1'  +  •  •  •  +  e^' v]. 


(2-11) 


Subtracting  Equation  (2-11)  from  Equation  (2-10)  gives 


or 


AF(1  -  ejv )  =  A0  (1  -  e  ^  v ) 

I 


AF  =  A0 


l-ejNeV 

l-ejv 


(2-12) 

(2-13) 


This  may  be  rewritten  in  a  more  convenient  form  as  follows 


e^-l  pjNeV/2f  jNey/2  -jNey/2*] 

AF  =  A  - _ _  =  a  - _  e  ~e 

0  ejv-l  0  ejv/2  eiv/2_e-jy/2  ’ 


so  that 
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(2-14) 


AF  =  A 


j(N.-l)V/2 

oe 


sin(Ne\j//2) 

sin(\|r/2) 


The  phase  factor  of  the  array  is  given  by  e)(N'~1)¥/2.  If  the  array  was  centered  about  the 
origin,  this  phase  factor  would  not  be  present  since  it  represents  the  phase  shift  of  the  array 
phase  center  relative  to  the  origin.  Neglecting  the  phase  factor,  Equation  (2-14)  becomes 


AF  =  A0 


sin(NeV}//2) 

sin(\|//2) 


(2-15) 


The  maximum  value  of  the  array  factor  is  equal  to  A0Ne  and  occurs  when  \|/=0.  To 
normalize  it  so  that  the  maximum  value  is  equal  to  unity  (if  A0=l)  and  is  independent  of  the 
number  of  elements,  Equation  (2-15)  is  divided  by  Ne 


1  sm(Ney/2) 
Ne  sin(\j//2) 


(2-16) 


Now  let  0O  be  the  corresponding  value  of  0  for  which  the  array  factor  is  maximum.  Then 
from  Equation  (2-7)  we  have  \\r  =  kd  sin0o  +  a  =  0 ,  or 


a  =  -kdsin90  = 


2n  d 
~ 


sin0o. 


(2-17) 


Equation  (2-17)  represents  the  element-to-element  phase  shift  in  the  excitation  currents 
necessary  to  produce  an  array  factor  main  beam  maximum  in  a  direction  60  relative  to  the 
line  along  which  the  array  elements  are  aligned.  This  main  beam  scanning  by  phase  control 
mode  can  be  implemented  by  substituting  Equation  (2-17)  into  Equation  (2-7)  obtaining 
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V  =  kd(sin  0  -  sin  90 )  =  (sin  0  -  sin  0O ) . 

A/ 


(2-18) 


Setting  the  phases  of  the  array  in  this  way  is  called  “steering  the  beam”  [Ref.  3]. 


C.  TWO-ELEMENT  INTERFEROMETER  ANTENNA 

For  a  two-element  interferometer  antenna  (Ne=2),  Equation  (2-16)  becomes 


AF  '  -  1  sinW 
nonn  2  sin(\|//2)' 


Assuming  the  two-element  configuration  is  centered  about  the  origin  of  the  reference 
system,  and  using  the  trigonometric  identity 


sin(2\|/)  =  2sin(v)cos(v), 


Equation  (2-19)  can  be  simplified  further  becoming 


AF„, 


sin(\j//2)cos(y/2) 

sin(\jr/2) 


or 


AFno™  =  cos(y/2). 

A  plot  of  the  array  factor  using  Equation  (2-21)  is  shown  in  Figure  2.4  for  d=6L 


(2-19) 


(2-20) 


(2-21) 
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Normalized  AF  Amplitude 


i 


Figure  2.4:  Array  Factor  Amplitude  for  a  Two-Element  Linear  Array  Separated  by  a 

Distance  d=6L 


The  phase  response  of  the  two-element  interferometer  antenna  with  d=6A,  is  derived  from 
Equation  (2-18)  with  0o=O  and  -90°  <  0  <  90°  as  illustrated  in  Figure  2.5a.  The  sawtooth 
waveform  results  by  constraining  \j/  between  -ic  and  n. 
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Phase  Response  ( rad ) 


Figure  2.5a:  Phase  Response  of  a  Two-Element  Linear  Array  Separated  by  a 

Distance  d=6  X. 


Figure  2.5b  shows  the  resultant  phase  response  folding  waveform  that  will  be  encoded 
using  the  RNS. 
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Figure  2.5b:  Phase  Response  of  a  Two-Element  Linear  Array  Separated  by  a  Distance 
d=6  X  (\|/  constrained  to  vary  between  -7t  and  7t). 


Note  the  sawtooth  nature  of  this  folding  waveform.  Also  note  the  change  in  folding  period 
as  the  input  direction  of  arrival  approaches  ±90°. 
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III.  SYMMETRICAL  NUMBER  SYSTEM  ENCODING  OF  THE 

AMPLITUDE  RESPONSE 


A.  THE  SYMMETRICAL  NUMBER  SYSTEM  (SNS) 

The  optimum  SNS  is  composed  of  a  number  of  Pairwise  Relatively  Prime  (PRP) 
moduli  nij.  The  integers  within  each  SNS  modulus  are  derived  by  amplitude  analyzing  a 
symmetrically  folded  waveform  with  the  period  of  the  waveform  equal  to  twice  the  PRP 
modulus,  i.e.,  2vc\  [Ref.  5].  Given  m,  the  integer  values  within  twice  the  individual 
modulus  are  given  by  the  row  vector 


x  =  [0,l,-,m-l,m-l,-,l,0].  (3-1) 

From  (3-1)  the  required  number  of  threshold  levels  is  m;-l.  By  analyzing  Equation  (3-1)  it 
is  clear  that  the  integers  within  each  modulus  do  not  form  a  complete  system  of  length  2mi 
by  themselves  due  to  the  presence  of  ambiguities.  However,  the  ambiguities  that  arise 
within  the  modulus  can  be  resolved  by  using  an  N  moduli  set  (mp  m2,...,mN)  and 
combining  the  values  from  all  the  channels  together.  By  doing  this,  the  SNS  becomes  a 
complete  system  having  a  one-to-one  correspondence  with  the  residue  number  system.  For 
N  equal  to  the  number  of  PRP  moduli,  the  dynamic  range  of  this  scheme  is  [Ref.  5, 6] 

M  =  f[m,  (3-2) 

i=l 

This  dynamic  range  is  also  the  position  of  the  first  repetitive  moduli  vector  and  the  total 
number  of  quantization  levels  without  ambiguities.  As  an  example  consider  the  moduli  set 
mj=3, 1^=4,  and  m3=5.  The  first  repetitive  moduli  vector  occurs  at  an  input  equal  to  60, 
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and  the  number  of  quantization  levels  without  ambiguities  is  also  60  (0  to  59).  Table  3.1 
illustrates  the  SNS  for  the  moduli  set  (3,  4,  5).  Figure  3.1  shows  the  symmetrical  folding 
waveforms  and  SNS  states  for  mj=3,  m2=4,  and  m3=5.  The  individual  states  within  the 
SNS  represent  the  number  of  threshold  levels  crossed  by  the  folding  waveform.  A  total  of 
nyl  threshold  levels  are  required  for  each  folding  waveform. 


Normalized 

Input 

mj=3, 

Moduli 

m2=4, 

m3=5 

6  i 

0 

0 

0 

l 

1 

1 

1 

2 

2 

2 

2 

3 

2 

3 

3  • 

4 

1 

3 

4 

5 

0 

2 

4 

6 

0 

1 

3 

7 

1 

0 

2 

8 

2 

0 

1 

9 

2 

1 

0 

10 

1 

2 

0 

11 

0 

3 

1 

12 

0 

3 

2 

« 

52 

1 

3 

2 

53 

0 

2 

3 

54 

0 

1 

4 

55 

1 

0 

4 

56 

2 

0 

3 

57 

2 

1 

2 

58 

1 

2 

1 

59 

0 

3 

0 

60 

0 

3 

0 

Table  3.1:  An  Example  to  Illustrate  the  Optimum  SNS 
Using  Moduli  Set  1^=3, 1^=4,  m3=5. 
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Array  Factor  Output 


□  :  SNS  states  (number  of  threshold  levels  crossed) 


0  1  2  3  4  5  6  7  8  9  10  11  12  13  14  15 

Normalized  Input  Direction  of  Arrival 


Figure  3.1:  Optimum  SNS  Waveforms  and  Output  Codes  for  mt=3, 1x^=4,  and  m3=5. 


B.  SNS  ANTENNA  ARCHITECTURE 


Figure  3.2  shows  a  schematic  representation  of  the  SNS  antenna  architecture  for 
instantaneous  direction  finding.  The  array  consists  of  four  elements.  Element  1  is  combined 
with  elements  2,  3,  and  4.  The  three  pairs  of  elements  (1,  4),  (1,  3),  and  (1,  2)  form  a 
three-baseline  modular  interferometer  with  mj=3, 1x12=4  and  1x13=5  respectively. 


Figure  3.2:  Schematic  Representation  of  the  Proposed  SNS  Antenna 
Architecture  for  Combat  Direction  Finding. 
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The  array  factor  amplitude  response  from  each  pair  of  elements  (corresponding  to  modulus 
n^)  is  quantized  using  rr^-l  comparators.  The  comparator  threshold  levels  are  adjusted  to 
midlevel  quantize  the  folded  amplitude  response  into  the  SNS  format.  The  SNS-to-decimal 
logic  block  then  recombines  each  channel  and  converts  the  SNS  representation  to  a  more 
familiar  digital  output 

The  distance  between  each  pair  of  elements  must  be  determined  in  order  to  obtain 
the  correct  number  of  folds  needed  to  cover  the  entire  field  of  view  (from  0  to  ±90°)  for  the 
chosen  moduli  set.  That  is,  the  distance  between  each  pair  of  elements  must  be  derived  as  a 
function  of  the  PRP  modulus,  m;,  so  that  the  corresponding  amplitude  response  has  the 
correct  number  of  folds.  For  a  given  modulus  mi?  the  number  of  folds  needed  to  cover  the 
dynamic  range  is  given  by 


ni,SNS 


dynamic  range  M 
2  m;  2m;  ’ 


(3-3) 


where  M  is  given  by  Equation  (3-2).  Next  consider  the  array  factor  given  by  Equation  (2- 
21).  Define  E=\|//2  given  by 


£  =  —  sin9.  (3-4) 

From  Equation  (3-4)  we  have  that: 


e  =  o° 

=> 

e  =  9o° 

=> 

&  7C  d 

4=  l 

(3-5) 

CD 

II 

vb 

o 

=> 

y  K  d 

X  * 
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Each  fold  of  the  array  factor  corresponds  to  a  change  in  phase  of  2%  which  implies 

^(90° )  -  ^(-90°  )  =  2%.  (3-6) 

Substituting  in  Equation  (3-6)  the  values  of  £(±90°)  found  in  Equation  (3-5)  we  obtain 

2jtd  .  ,  . 

—  =  2%,  d  =  X.  (3-7) 

Therefore,  the  distance  between  elements  is  given  by 

M 
2m; 

where  X  is  the  wavelength  corresponding  to  the  center  frequency  used  to  construct  the 
antenna.  For  example,  if  an  antenna  is  to  be  designed  assuming  a  center  frequency  f0=8.5 
GHz  using  rn,=3, 1^=4  and  1^=5,  the  distance  required  between  each  pair  of  elements  is: 

dt  =  10  X  =  35.29  cm  (rrij  =  3) 
d2  =  7.5  X  =  26.47  cm  (m2  =  4) 
d3  =  6  A,  =21.18  cm  (m3  =  5). 

With  the  correct  distances  between  elements,  the  amplitude  response  from  each  pair  of 
elements  will  be  quantized  using  m^l  comparators.  The  function  of  the  comparators  is  to 
generate  the  integers  within  each  modulus  by  quantizing  the  amplitude  response  according 
to  the  modulus  n^  with  which  it  is  associated.  Finally  the  recombination  of  the  SNS  code 
from  each  pair  of  elements  will  be  the  instantaneous  direction  of  arrival. 
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C.  SIMULATION  RESULTS 


The  array  factor  plot  for  the  nij=3  array  as  the  input  direction  of  arrival  varies  from 
-90°  to  90°  is  shown  in  Figure  3.3.  As  shown  the  folding  period  changes  as  the  angle  6 
goes  from  0°  to  ±  90°. 


Figure  3.3:  Array  Factor  as  a  Function  of  0  for  a  Two-Element  Array  Antenna 

With  d=7.5  X. 


The  SNS  code  cannot  be  applied  to  this  waveform  in  a  straight  forward  manner.  However, 
if  the  array  factor  is  examined  as  a  function  of  sin(0)  as  shown  in  Figure  3.4,  the  array 
factor  corresponds  to  a  symmetrically  folded  sinusoidal  waveform  (-1  <  sin(0)  <  1)  with  a 
constant  period,  thereby  satisfying  the  principal  condition  for  applying  the  SNS  code. 

The  chosen  moduli  have  a  dynamic  range  M=60,  therefore  the  quantization  width, 
considering  the  field  of  view  (from  0°  to  ±  90°)  that  corresponds  to  the  variation  in  sin0 
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(from  0  to  ±  1),  will  be  equal  to  2/60=0.0333.  For  a  true  symmetrically  folded  sinusoidal 
waveform  as  a  function  of  0,  with  constant  period,  the  quantization  width  would  be 
constant,  independent  of  the  direction  of  arrival,  and  equal  to  3°  (180760).  This,  however, 
is  not  the  case.  Comparing  Figures  3.3  and  3.4  it  is  possible  never  the  less  to  establish  a 
relationship  between  the  constant  quantization  width  (when  the  array  factor  is  plotted  as  a 
function  of  sin(0))  and  the  corresponding  quantization  width  as  a  function  of  0  (when  the 
array  factor  is  plotted  as  a  function  of  0). 


Figure  3.4:  Array  Factor  as  a  Function  of  sin(0)  for  a  Two-Element  Array  Antenna  With 

d=7.5  l. 


As  a  function  of  0  the  quantization  width  depends  on  the  direction  of  arrival  of  the 
incoming  signal  and  is  not  uniform.  The  transition  points  from  one  quantization  level  to 
another  quantization  level  occur  at  an  angle  0t  given  by 
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0t  =  sin  :(x) 


(3-9) 


where  x  is  given  by 


x  =  -l,-l  +  |-,-l  +  -i  .,1.  (3-10) 

M  M 

Table  3.2  shows  the  resulting  transition  points,  0t,  and  the  respective  quantization  widths. 

The  array  factor  for  each  pair  of  elements  {(1, 4),  (1, 3),  and  (1, 2)}  corresponding 
to  modulus  mj=3,  m2=4,  and  m3=5  respectively  are  shown  in  Figures  3.5,  3.6,  and  3.7. 
Also  shown  are  the  states  within  each  modulus.  The  array  factor  is  phased  at  -90°  in  all 
cases  so  that  the  SNS  code  ‘000’  will  correspond  to  our  first  quantization  level  at  -90°  (see 
Table  3.1).  To  obtain  this  proper  phase  relationship  between  pairs,  a  phase  shift  equal  to  7t , 
Jt/2,  and  7t  was  added  to  elements  (1, 2),  (1, 3),  and  (1, 4)  respectively. 

The  next  step  is  to  separately  quantize  the  array  factor  waveform  corresponding  to 
each  pair  of  elements  using  the  SNS  encoding  procedure  explained  in  Section  A  of  this 
chapter.  The  number  of  threshold  levels  (or  comparators)  depends  on  the  modulus  rr^  that 
is  being  considered  and  is  given  by  m;-l  as  previously  noted.  The  matching  threshold 
values  for  the  comparators  are  found  by  simulation.  The  transition  points  from  one 
quantization  level  to  another  occur  at  the  values  of  0t  shown  in  Table  3.2. 
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sin1  x 

et 

(degrees) 

Quantization 

Width 

sin'1  x 

et 

(degrees) 

Quantization 

Width 

-1.0000 

-90.0000 

0.0000 

0.0000 

1.9102 

-0.9666 

-75.1649 

14.8351 

0.0333 

1.9102 

1.9123 

-0.9333 

-68.9605 

6.2044 

0.0666 

3.8102 

1.9166 

BIKIIIIII1 

-64.1581 

4.8025 

0.1000 

5.7392 

1.9231 

-0.8666 

-60.0736 

4.0845 

0.1333 

7.6623 

1.9318 

-0.8333 

-56.4427 

3.6309 

0.1666 

9.5942 

1.9429 

-0.8000 

-53.1301 

3.3126 

0.2000 

11.5370 

1.9564 

-0.7666 

-50.0555 

3.0746 

0.2333 

13.4934 

1.9726 

-0.7333 

-47.1666 

2.8889 

0.2666 

15.4660 

1.9916 

-0.7000 

-44.4270 

2.7396 

0.3000 

17.4576 

2.0136 

-0.6666 

-41.8103 

2.6167 

0.3333 

19.4712 

2.0390 

-0.6333 

-39.2965 

2.5138 

0.3666 

21.5102 

2.0680 

-0.6000 

-36.8699 

2.4266 

0.4000 

23.5782 

2.1011  | 

-0.5666 

-34.5181 

2.3518 

0.4333 

25.6793 

2.1389 

-0.5333 

-32.2310 

2.2872 

0.4666 

27.8181 

2.1819 

-0.5000 

-30.0000 

2.2310 

0.5000 

30.0000 

2.2310  t 

-0.4666 

-27.8181 

2.1819 

0.5333 

32.2310 

2.2872 

-0.4333 

-25.6793 

2.1389 

0.5666 

34.5181 

2.3518 

-0.4000 

-23.5782 

2.1011 

0.6000 

36.8699 

2,4266 

-0.3666 

-21.5102 

2.0680 

0.6333 

39.2965 

2.5138 

-0.3333 

-19.4712 

2.0390 

0.6666 

41.8103 

2.6167 

-0.3000 

-17.4576 

2.0136 

0.7000 

44,4270 

2.7396 

-0.2666 

-15.4660 

1.9916 

0.7333 

47.1666 

2.8889 

-0.2333 

-13.4934 

1.9726 

0.7666 

50.0555 

3.0746 

-0.2000 

-11.5370 

1.9564 

0.8000 

53.1301 

3.3126 

-0.1666 

-9.5942 

1.9429 

0.8333 

56.4427 

3.6309 

-0.1333 

-7.6623 

1.9318 

0.8666 

60.0736 

4.0845 

-0.1000 

-5.7392 

1.9231 

0.9000 

64.1581 

4.8025 

-0.0666 

-3.8102 

1.9166 

0.9333 

68.9605 

6.2044 

-0.0333 

-1.9102 

1.9123 

0.9666 

75.1640 

14.8351 

Hi 

0.0000 

1.9102 

1.0000 

90.0000 

Table  3.2:  Equivalent  Quantization  Levels  Transition  Points  and  Quantization 

Level  Width. 
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Figure  3.5:  Array  Factor  Waveform  and  Respective  SNS  Quantization  Code  for  the  Pair  of 
Elements  (1,4)  Corresponding  to  Modulus  mj=3. 


Figure  3.6:  Array  Factor  Waveform  and  Respective  SNS  Quantization  Code  for  the  Pair  of 
Elements  (1,3)  Corresponding  to  Modulus  m2=4. 
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Figure  3.7 :  Array  Factor  Waveform  and  Respective  SNS  Quantization  Code  for  the  Pair  of 
Elements  (1,2)  Corresponding  to  Modulus  m3=5. 


The  normalized  comparator  matching  threshold  values  for  each  pair  of  elements 
corresponding  to  modulus  3, 4,  and  5  are  shown  in  Table  3.3. 


Elements 

Modulus 

T1 

T2 

T3 

T4 

(1,2) 

5 

-0.8090 

-0.3090 

0.3090 

0.8090 

(1,3) 

4 

-0.7071 

0.0000 

0.7071 

(1,4) 

3 

-0.5000 

0.5000 

i 

i 

i 

i 

i 

t 

Table  3.3:  Normalized  Comparator  Matching  Threshold  Values  for  Each  Pair  of  Elements. 
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The  comparator  states  (or  thermometer  codes)  from  each  channel  are  recombined  in 
an  SNS  logic  block  to  give  the  instantaneous  resolved  direction  of  arrival.  Figure  3.8 
shows  the  resolved  direction  of  arrival  or  antenna  transfer  function  corresponding  to  the  60 
different  possible  combinations  of  the  SNS  for  the  moduli  set  (3,  4,  5).  Note  that  the 
quantization  width  is  not  constant 


Figure  3.8:  SNS  Resolved  Direction  of  Arrival  for  -90°  <  0  <  90°(no  element  factor). 


The  actual  amplitude  response  of  a  linear  array  antenna  is  given  by  Equation  (2-1) 
which  means  that  a  multiplicative  factor  (the  element  factor)  has  to  be  taken  into 
consideration  in  addition  to  the  array  factor.  The  effect  of  the  element  factor  is  to  introduce 
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a  drop-off  in  the  antenna  gain.  The  element  factor  will  depend  on  the  type  of  elements  that 
are  to  be  used  in  the  array  antenna  (waveguides,  dipoles,  etc.).  For  example  consider  that 
this  element  factor  is  characterized  by  a  cosine  waveform.  Taking  this  characterization  into 
account.  Equation  (2-1)  can  be  rewritten  as  follows 

yjcos(e).  (3-11) 


E(0, 4>  =0)  =  A0  cos 


Figure  3.9  shows  a  plot  of  the  electrical  field  response  of  the  array  assuming  that  the 
element  factor  is  characterized  by  cos(0).  The  element  factor  modulates  the  array  factor  so 
that  the  resultant  waveform  no  longer  has  a  constant  amplitude. 


Figure  3.9:  Amplitude  of  the  Electrical  Field  Response  of  a  Two  Linear  Array:  Elements  1 

and  2. 
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When  the  SNS  quantization  methodology  is  applied  to  the  resultant  waveform  (same 
threshold  levels)  a  number  of  errors  result  as  shown  in  Figure  3.10.  However,  the  results 
are  still  acceptable  in  a  smaller  interval  (-20°  <  0  <  20°).  The  antenna  transfer  function  for 
this  narrow  region  is  shown  in  Figure  3.11.  The  encoding  errors  which  exist  at  the  code 
transition  points  are  the  result  of  some  of  the  comparators  not  switching  simultaneously  as 
they  should.  These  errors  are  easily  isolated  by  using  a  few  additional  comparators  in  the 
minimum  modulus  channel  [Ref.  5].  Investigation  into  the  error  isolation  however,  is 
beyond  the  scope  of  this  thesis  and  is  subject  of  future  research. 


Figure  3. 10:  SNS  Resolved  Direction  of  Arrival  for  -90°  <  0  <  90°  (element  factor 

included). 


Figure  3.1 1:  SNS  Resolved  Direction  of  Arrival  for  -20°  <  0  <  20°(element  factor 

included). 


D.  EXPERIMENTAL  RESULTS 

A  linear  antenna  was  built  based  on  the  theory  developed  previously  using  four 
rectangular  waveguides  mounted  in  a  large  ground  plane.  Figure  3.12  shows  the 
configuration  of  the  antenna  (front  plane,  back  plane,  and  mounted  on  the  pedestal  in  the 
anechoic  chamber).  The  connection  between  elements  are  as  illustrated  in  Figure  3.2. 
Using  the  anechoic  chamber  facility  at  the  Naval  Postgraduate  School,  antenna  pattern 
measurements  (power  response)  for  the  H-field  were  taken  for  each  pair  of  elements. 
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Figure  3.12.b:  Photograph  of  the  Built  Antenna  (Back  Plane). 
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Figure  3.12.c:  Photograph  of  the  Built  Antenna  Mounted  on  the  Pedestal  in  the 

Anechoic  Chamber. 


Before  evaluating  the  experimental  results,  the  simulated  results  for  this  particular  antenna 
(rectangular  waveguide  elements)  are  examined.  The  theoretical  element  factor  (EF)  for  a 
rectangular  waveguide  when  analyzing  the  H-field  is  given  by 


EF  -  cos(ka/2  sin  6) 
it2  -(ka/2  sin0)2 


COS0 


(3-12) 


where  a  is  the  width  of  the  waveguide  and  k  is  the  wave  number  [Ref.  3].  When  Equation 
(3-12)  is  used  as  the  element  factor  (instead  of  the  cos(0)  used  previously  in  Section  C), 
the  resultant  electrical  field  for  each  pair  of  elements  is  shown  in  Figure  3.13.  Figure  3.14 
shows  the  resolved  direction  of  arrival  obtained  when  this  electrical  field  is  processed  using 
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Electrical  Field  Response 


the  SNS  architecture  presented  in  Section  A  of  this  chapter.  The  results  compare  favorably 
to  those  shown  in  Figure  3-10. 


Input  Direction  of  Arrival 

Figure  3.13:  Electrical  Field  Response  for  Mod  3,  Mod  4,  and  Mod  5  When  Using 

Rectangular  Waveguide  Elements. 


37 


Figure  3.14:  SNS  Resolved  Direction  of  Arrival  as  a  Function  of  0,  When  Using  the  EF 

for  a  Rectangular  Waveguide. 


The  experimental  measurements  of  the  antenna  pattern  for  each  pair  of  elements  is 
shown  in  Figure  3.15.  Since  the  SNS  architecture  quantizes  the  amplitude  of  the  received 
electrical  field  (not  the  power)  it  is  necessary  to  convert  to  the  collected  measurements 
corresponding  electrical  field  response.  The  relationship  between  the  Power  (P)  and  the 
Electrical  Field  (E)  is  given  by 


P  =  20  Log10 1 E 1. 


(3-13) 
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Antenna  Pattern  (  dB  ) 


The  absolute  value  of  the  electrical  field  is  then 


I E  l  =  10p/2°. 


(3.14) 


Input  Angle  of  Arrival 


Figure  3.15:  Antenna  Pattern  for  Each  Pair  of  Elements  (Experimental  Data). 


Using  Equation  (3-14),  the  experimental  measurements  for  each  modulus  (pair  of  elements) 
were  used  to  calculate  the  absolute  value  of  the  electrical  fields  and  are  shown  in  Figure 
3.16. 
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Absolute  value  of  the  Electrical  Field  Response  Normalized 


Electrical  Field  Response  ( volts) 


1 


Figure  3.17.  Electrical  field  Response  for  Each  Pair  of  Elements  Corresponding  to 

Modulus  3,  4,  and  5. 


The  electrical  field  measurements  for  each  pair  of  elements,  representing  modulus  3,  4,  and 
5  is  then  quantized,  and  the  output  from  each  channel  is  recombined  to  obtain  the  resolved 
direction  of  arrival.  The  resolved  direction  of  arrival  is  shown  in  Figure  3.18. 
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Figure  3.18:  Resolved  Direction  of  Arrival  Using  Experimental  Data. 


Comparing  the  experimental  results  (Figure  3.18)  with  the  simulated  results  (Figure  3.14) 
shows  a  good  deal  of  agreement.  However,  considering  just  the  interval  from  0°  to  ±  20°, 
we  notice  many  more  encoding  errors  in  the  experimental  results.  If  the  electrical  field  was 
available  directly  in  the  anechoic  chamber,  the  approximation  would  not  have  been 
necessary  and  the  number  of  encoding  errors  would  have  been  more  in  line  with  those  of 
Figure  3.14. 
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IV.  RESIDUE  NUMBER  SYSTEM  ENCODING  OF  THE  PHASE 

RESPONSE 

A.  THE  RESIDUE  NUMBER  SYSTEM  (RNS) 

The  ancient  study  of  the  Residue  Number  System  begins  with  a  verse  from  a  third- 
century  book,  Suan-ching,  by  Sun  Tzu  [Ref.  7]: 

We  have  things  of  which  we  do  not  know  the  number. 

If  we  count  them  by  threes,  the  remainder  is  2, 

If  we  count  them  by  fives,  the  remainder  is  3, 

If  we  count  them  by  sevens,  the  remainder  is  2, 

How  many  things  are  there? 

The  answer,  23. 

This  1700-year-old  number  system  has  been  attracting  a  great  deal  of  attention 
recently.  Digital  systems  that  use  residue  arithmetic  units  may  play  an  important  role  in 
ultra-high  speed,  dedicated,  real-time  systems  that  support  pure  parallel  processing  of 
integer  valued  data.  It  is  a  “carry  free”  system  that  performs  addition,  subtraction,  and 
multiplication  as  concurrent  operations,  side-stepping  one  of  the  principal  arithmetic  delays: 
managing  carry  information.  Some  applications  such  as  fast  Fourier  transform,  digital 
filtering,  and  image  processing  utilize  the  efficiencies  of  RNS  arithmetic  in  addition  and 
multiplication.  The  RNS  has  been  employed  efficiently  in  the  implementation  of  digital 
processors. 

Similar  to  the  SNS,  the  RNS  is  also  composed  of  a  number  of  pairwise  relatively 
prime  moduli  rrij  [Ref.  6].  The  integers  within  each  RNS  modulus  are  representative  of  a 
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saw-tooth  waveform  with  the  period  of  the  waveform  equal  to  the  modulus.  In  this 
scheme,  given  m,  the  integer  values  within  the  individual  modulus  are  given  by  the  row 
vector 


y  =  [0,l,...,m-l].  (4-1) 

Again,  the  required  number  of  comparators  for  each  channel  is  m-l.  Notice  that  all  the 
integers  within  Equation  (4-1)  form  a  complete  system  of  length  (no  ambiguities).  For  N 
equal  to  the  number  of  moduli,  the  dynamic  range  of  this  scheme,  as  well  as  the  number  of 
quantization  levels,  is  given  by 


M  =  fK  (4-2) 

i=l 


Table  4.1  illustrates  the  RNS  for  the  moduli  1^=3,  m2=4,  m3=5.  Figure  4.1  shows  the 
RNS  folding  waveforms  and  the  corresponding  RNS  states. 


44 


7  13  2 

8  2  0  3 

9  0  14 

10  1  2  0 

11  2  3  1 

12  0  0  2 


52  1 

53  2 

54  0 

55  1 

56  2 

57  0 

58  1 

59  _ 2 

60  0 


Table  4.1:  An  Example  to  Illustrate  the  RNS  Encoding  Procedure  for 
the  Moduli  1^=3, 1112=4  and  m3=5. 
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□  :  RNS  states  (threshold  levels  crossed) 


0  1  2345678  9  10  11 

Normalized  Input  Direction  of  Arrival 

Figure  4.1:  RNS  Folding  Waveforms  and  Output  Codes  for  the  Moduli 

0^=3, 1112=4,  and  m3=5. 

B.  RNS  ANTENNA  ARCHITECTURE 

The  phase  response  of  a  two-element  interferometer  corresponds  to  a  saw-tooth 
folding  waveform.  The  RNS  antenna  architecture  being  investigated  is  shown  in  Figure 
4.2.  Phase  detectors  are  used  to  combine  the  individual  pairs  of  elements.  Also,  the 
comparator  threshold  levels  and  logic  block  are  particular  to  the  RNS  (see  Table  4.1). 
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DIRECTION  OF  ARRIVAL  0 


Figure  4.2:  RNS  Antenna  Architecture. 


The  distance  between  each  pair  of  elements  must  be  derived  as  a  function  of  the 
PRP  modulus,  nr,  so  that  the  corresponding  phase  response  has  the  correct  number  of 
folds.  For  a  given  modulus  m*,  the  number  of  folds  needed  to  apply  the  RNS  code  within 
the  dynamic  range  is  given  by 
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ni,RNS  " 


_  dynamic  range  _  M 


m- 


nr 


(4-3) 


Next  consider  the  two  element  interferometer  shown  in  Figure  1.1.  The  phase  difference 
between  the  received  signal  due  to  the  separation  between  elements  is  given  by 


A<J)  = 

27td  . 

=  — —  sin0. 

X 

(4-4) 

From  Equation  (4-4)  we  have  that: 

e 

=  0° 

=>  A(J)  =  0 

e 

=  90° 

A ,  2nd 

(4-5) 

e 

=  -90° 

-  A*=-r. 

Since  each  fold  of  the  phase  response  corresponds  to  a  change  in  phase  of  2jt,  we  require 
that 


A<|)(90°)-A<|>(-90°)  =  27t  (4-6) 

which  implies  (from  Equation  (4-5)  , 


Finally,  the  distance  between  elements  is  given  by 
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(4-8) 


M 

mi 


Using  the  same  moduli  set  with  center  frequency  of  f0=8.5  GHz,  the  distances  for  each  pair 
of  elements  are 


dj  =10  X, 
d2  =  7.5  X, 
d3  =  6.0  X. 

which  are  the  same  as  those  given  for  the  SNS  architecture.  That  is,  the  same  hardware 
constructed  for  the  SNS  array  can  be  used  for  the  RNS  array. 

C.  SIMULATION  RESULTS 

The  correct  distances  between  each  pair  of  elements  were  found  in  Section  B  so  that 
the  output  waveforms  representative  of  modulus  0^=3,  m2=4,  and  m3=5,  have  the  correct 
number  of  folds  required  to  quantize  the  entire  phase  response  within  the  field  of  view 
(from  0°  to  ±90°).  Figure  4.3  shows  a  plot  of  the  phase  response,  as  a  function  of  the  input 
direction  of  arrival  for  a  two-element  interferometer  separated  by  a  distance  d=6A,.  The 
saw-tooth  waveform  stretches  out  as  0,  the  input  direction  of  arrival,  goes  from  0°  to  ±90°. 
However,  if  the  phase  response  is  plotted  as  a  function  of  sin(0),  a  perfect  saw-tooth 
waveform  is  obtained  as  shown  in  Figure  4.4. 
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Phase  Response  (rad)  Phase  Response  (rad) 


Figure  4.3:  Phase  Response  as  a  Function  of  the  Input  Direction  of  Arrival  for  a 
Two-Element  Interferometer  with  Separation  d3=6X,. 


Figure  4.4:  Phase  Response  as  a  Function  of  sin(0)  for  a  Two-Element  Interferometer  with 

Separation  d3=6X,. 
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Based  on  these  results,  the  same  analysis  that  was  performed  when  dealing  with  the  array 
factor  can  be  applied  here.  In  Figures  4.5,  4.6,  and  4.7,  the  phase  response  waveform  and 
respective  output  quantization  code  for  modulus  mj=3,  modulus  m^,  and  modulus  m3=5 
are  plotted,  as  a  function  of  the  input  direction  of  arrival. 


Figure  4.5:  Phase  Response  Waveform  and  Respective  RNS  Quantization  Code  for  the 
Pair  of  Elements  (1,4)  Corresponding  to  Modulus  mj=3. 
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Figure  4.6:  Phase  Response  Waveform  and  Respective  RNS  Quantization  Code  for  the 
Pair  of  Elements  (1,3)  Corresponding  to  Modulus  0^=4. 


Figure  4.7:  Phase  Response  Waveform  and  Respective  RNS  Quantization  Code  for  the 
Pair  of  Elements  (1,2)  Corresponding  to  Modulus  1^=5. 
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The  quantization  of  the  phase  response  from  each  pair  of  elements  is  made 
separately  using  a  small  bank  of  comparators.  The  number  of  threshold  levels  (and 
comparators)  are  the  same  as  that  for  the  array  factor  amplitude  analysis.  However,  the 
matching  threshold  values  are  different  as  shown  in  Table  4.2. 


Elements 

Modulus 

T1 

T2 

T3 

T4 

(1,2) 

5 

-1.8850 

-0.6283 

0.6283 

1.8850 

(1,3) 

4 

-1.5708 

0.0000 

1.5708 

(1,4) 

3 

-1.0472 

1.0472 

Table  4.2:  Comparator  Threshold  Levels  for  Each  Pair  of  Elements  When 
Considering  the  Phase  Analysis. 


For  a  specific  input  direction  of  arrival  and  pair  of  elements,  there  will  be  a 
corresponding  integer  (number  of  comparators  in  the  on  state)  within  the  PRP  modulus  m,. 
The  value  of  those  integers  depends  on  the  measured  phase  difference  associated  with  each 
pair  of  elements  for  that  specific  input  direction  of  arrival.  The  comparator  state  for  each 
pair  together  can  then  be  represented  by  3  RNS  codes.  The  3  RNS  codes  are  indexed  to  an 
RNS  logic  block  to  give  the  resolved  direction  of  arrival.  Table  4.3  shows  the  3  RNS 
codes  and  the  resolved  direction  of  arrival. 
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3  RNS 
Codes 


000 


111 


222 


033 


104 


210 


021 


132 


203 


014 


120 


231 


002 


113 


224 


030 


101 


212 


023 


134 


200 


Oil 


122 


233 


004 


110 


221 


032 


103 


214 


Resolved  DOA 
(degrees) 


10.5655 


31.1155 


33.3745 


40.5534 


45.7968 


62.1158 


3  RNS 

Resolved  DOA 

Code 

(degrees) 

020 

131 

1  -72.0627  1 

202 

[ 

013 

124 

-58.2581  1 

230 

—  Lilli  1 LH 

001 

1  -51.5928  | 

112 

■BSsUlSH 

223 

WKBSSSsMM 

034 

100 

i 

211 

-38.0832 

022 

-35.6940 

133 

204 

■KflDEm 

010 

121 

WKKSSsSSSSUKM 

232 

-24.6287  1 

003 

114 

220 

-18.4644 

031 

-16.4618 

102 

-14.4797 

213 

024 

130 

201 

012 

123 

-2.8664  f 

234 

Table  4.3:  RNS  Code  and  Equivalent  Resolved  Direction  of  Arrival  for 
Modulus  3, 4,  and  5  Phased  at  0°  with  Center  Frequency  f0=8.5  GHz. 
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Figure  4.8,  shows  a  plot  of  the  resolved  direction  of  arrival  as  a  function  of  0,  the 
input  direction  of  arrival.  This  new  system  gives  the  direction  of  arrival  information  over 
the  complete  field  of  view  without  ambiguity.  However  the  accuracy  varies  depending  on 
the  actual  direction  of  arrival  of  the  incoming  signal.  Notice  also  that  the  quantization  width 
increases  as  the  direction  of  arrival  approaches  ±90°.  The  quantization  error  will  be 
addressed  in  Chapter  V. 


- 1 - 1 - 1 - 1 - 1 - 1 _ : _ I _ 1 _ L_ 

-80  -60  -40  -20  0  20  40  60  80 

Input  Direction  of  Arrival  (  degrees  ) 


Figure  4.8:  Resolved  Direction  of  Arrival  in  Function  of  the  Input  Direction  of  Arrival. 
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D.  EXPERIMENTAL  RESULTS 


Unfortunately  it  was  not  possible  to  obtain  acceptable  experimental  results  for  the 
RNS  antenna  architecture.  The  main  problem  was  the  nonavailability  of  the  phase 
detectors.  Another  problem  was  matching  the  elements  with  the  correct  length  of  cable  in 
order  to  obtain  the  correct  phase  difference  between  each  pair  of  elements.  Future 
investigations  are  currently  underway  and  will  correct  these  problem. 
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V.  FURTHER  CONSIDERATIONS 


A.  OVERVIEW 

Comparing  the  simulated  results  for  the  SNS  antenna  design  with  the  results 
obtained  for  the  RNS  antenna  design,  it  is  clear  that  the  second  case  permits  acquiring  the 
direction  of  arrival  information  over  the  complete  field  of  view  from  0°  to  ±90°.  When 
using  the  SNS  antenna  it  is  impossible  to  obtain  the  direction  of  arrival  information  over  the 
entire  range  from  -90°  to  90°  due  to  the  degradation  in  the  amplitude  response.  Since  the 
phase  response  has  a  “uniform”  distribution  (i.e.,  maximum  and  minimum  with  a  constant 
value)  over  the  complete  field  of  view,  it  is  possible  to  use  the  RNS  code  to  quantize  that 
signal  correctly.  Accordingly,  only  the  phase  analysis  will  be  considered  hereafter  since  the 
results  shown  are  more  encouraging.  At  this  point  we  have  not  presented  any  results 
concerning  the  following  aspects: 

•  quantization  errors  associated  with  this  new  antenna  design, 

•  the  antenna  frequency  response  (incoming  signal  frequency  different  than  the  one 
used  in  the  design), 


•  the  required  size  of  the  antenna  related  to  the  moduli  set  and  the  frequency  used  to 
build  die  antenna,  and 


•  the  influence  of  the  chosen  moduli  set  on  the  accuracy  and  quantization  error  of  the 
antenna. 

These  aspects  are  considered  in  this  chapter  for  the  RNS  architecture  (using  the  phase 
response).  Some  antennas  design  examples  are  given  to  show  the  degree  of  accuracy  that 
can  be  reached  with  this  direction  finding  architecture. 
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B.  QUANTIZATION  ERROR  ASSOCIATED  WITH  THE  RNS 
ARCHITECTURE 

The  conversion  of  an  analog  (continuous)  signal  into  a  digital  (discrete) 
representation  is  called  quantization.  Graphically,  the  quantization  process  means  that  a 
straight  line  representing  the  relation  between  the  input  and  output  of  a  linear  continuous 
system  is  replaced  by  a  staircase  characteristic  [Ref.  8],  as  shown  in  Figure  4.8.  The 
difference  of  two  adjacent  discrete  values  is  called  the  quantum  or  step  size.  Signals  applied 
to  a  quantizer,  with  the  input-output  characteristic  of  Figure  4.8,  are  sorted  into  amplitude 
slices,  and  all  input  signals  within  plus  or  minus  half  a  quantum  step  of  the  midvalue  of  a 
slice  are  replaced  in  the  output  by  the  midvalue.  The  quantizer  error  consists  of  the 
difference  between  the  input  and  output  signals  of  the  quantizer  so  that  the  maximum 
instantaneous  value  of  this  error  is  half  of  one  quantum  step,  and  the  total  range  of  variation 
is  from  minus  half  a  step  to  plus  half  a  step.  By  analyzing  Figure  4.8  again,  it  is  clear  that 
the  quantization  error  for  the  RNS  antenna  is  not  constant  since  the  code  widths  are  not 
constant.  The  Appendix  contains  the  Matlab™  program  code  to  calculate  the  quantization 
error  for  a  predetermined  moduli  set  as  a  function  of  the  input  direction  of  arrival.  Figure 
5.1  shows  a  plot  of  the  characteristic  quantization  error  curve  associated  with  the  resolved 
direction  of  arrival  from  Figure  4.8  for  the  moduli  set  (3, 4, 5)  and  for  a  frequency  equal  to 
8.5  GHz.  For  this  case,  the  results  show  that  the  quantization  error  varies  from  about  ±1° 
over  the  input  range  ±20°,  to  a  maximum  of  ±7.42°  for  direction  of  arrivals  near  ±90°. 
Figure  5.2  it  is  a  plot  of  both  the  resolved  direction  of  arrival  and  the  associated 
quantization  error  as  a  function  of  the  input  direction  of  arrival. 
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Figure  5.1:  Characteristic  Quantization  Error  Associated  With  the  RNS  Antenna  When 
Considering  the  Moduli  Set  (3, 4,  5)  and  the  f =f0=8.5  GHz. 


Figure  5.2:  Resolved  Direction  of  Arrival  and  Quantization  Error. 


C.  RNS  ANTENNA  FREQUENCY  RESPONSE 

A  frequency  (f0)  was  chosen  to  calculate  the  required  distance  between  each  pair  of 
elements  so  that  the  correct  number  of  folds  (for  each  pair  of  elements)  were  available  to 
apply  the  RNS  code  and  process  (quantize)  correctly  the  phase  response  without 
ambiguities  over  the  entire  field  of  view.  The  simulated  results  shown  previously  are 
obtained  assuming  that  the  frequency  of  the  incoming  signal  (fr)  was  equal  to  f0.  In  this 
case,  the  only  error  present  is  the  quantization  error.  It  is  interesting  to  examine  what 
happens  when  the  frequency  of  the  incoming  signal  has  some  other  value  different  from  f0. 
Figure  5.3,  shows  the  resolved  direction  of  arrival  for  an  incoming  signal  with  a  frequency 
f=4.25  GHz.  Comparing  these  results  with  the  results  in  Figure  4.8,  we  notice  the 
following:  the  total  number  of  quantization  levels  is  reduced  by  half  and  the  resolved 
direction  of  arrival  has  large  errors  as  the  input  direction  of  arrived  goes  from  0°  to  ±90°. 


Figure  5.3:  Resolved  Direction  of  Arrival  for  an  Incoming  Signal  With  a  Frequency 

f  =4.25  GHz. 
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Alternatively,  Figure  5.4  shows  the  resolved  direction  of  arrival  when  the  incoming  signal 
frequency  fr=17  GHz.  In  this  case,  we  notice  that  the  total  number  of  quantization  levels 
doubles  when  compared  with  the  simulated  results  shown  in  Figure  4.8  (fr=f0).  As  a  result, 
the  resolved  directions  of  arrival  are  ambiguous  because  the  same  3  RNS  codes  correspond 
to  more  that  one  direction  of  arrival.  In  addition  the  overall  error  is  very  high. 


Figure  5.4:  Resolved  Direction  of  Arrival  for  an  Incoming  Signal  With  a  Frequency 

fr  =  17  GHz. 


The  RNS  antenna  was  designed  assuming  a  moduli  set  (3, 4, 5)  and  a  frequency  (f0)  equal 
to  8.5  GHz.  For  a  given  PRP  modulus  m^  Equation  (4-3)  gives  the  required  number  of 
folds  to  apply  the  RNS  code  within  the  dynamic  range  without  ambiguities.  It  was  found  to 
be  equal  to  20  for  modulus  3.  In  addition,  the  distance  between  each  pair  of  elements  is 
given  by  Equation  (4-4).  This  distance  is  fixed  when  the  antenna  is  built  On  other  hand,  in 
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the  presence  of  a  signal  with  a  frequency  other  than  f0,  and  from  Equation  (1-1),  the  phase 
difference  between  the  received  signal  due  to  the  separation  between  elements  is  given  by 

A  ())  =  —  sin  9.  (5-1) 

Consider  modulus  3,  d=10A,0.  Substituting  that  value  in  Equation  (5-1)  gives 

A4,=  20£Alsin9.  (5-2) 

To  analyze  Equation  (5-2),  three  distinct  cases  need  to  be  considered  according  to  the  value 
of  the  ratio  X0/X[. 

•  Case  1:  XJXt  =  1  =>  fr  =  fr 

This  can  be  considered  the  ideal  case  since  the  corresponding  phase  response  (sawtooth 
waveform)  will  have  the  exact  number  of  folds  required  to  apply  the  RNS  code  within  the 
dynamic  range  without  ambiguities  over  the  complete  field  of  view  (from  0°  to  ±90°). 

•  Case  2:  XJXr  <  1  =>  fr<f#. 

In  this  case  the  phase  response  waveform  will  have  less  folds  than  the  required  number  to 
apply  the  RNS  code  within  the  dynamic  range.  No  ambiguities  will  be  present;  however, 
the  quantization  error  will  increase  since  a  smaller  number  of  quantization  levels  will  exist 
to  cover  the  same  field  of  view  (from  0°  to  ±90°).  This  case  is  directly  related  to  Figure  5.3 
where  the  ratio  XJX^.5.  The  equivalent  number  of  folds  and  quantization  levels  are  one 
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Phase  Response 


half  those  found  for  the  ideal  case.  Figure  5.5  shows  the  phase  response  for  both 
situations:  f=8.5  GHz  and  fr=4.25  GHz  for  modulus  m3=5. 
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Input  Direction  of  Arrival  (  degrees  ) 

Figure  5.5:  Phase  Response  Waveform  (Mod  5)  for  f  =8.5  GHz  (Dashed)  and 

f  =4.25  GHz. 

•  Case  3:  XJXr  >1  =>  fr  >  f#. 

For  this  case  the  resultant  phase  response  will  have  more  folds  than  the  required  number  to 
apply  the  RNS  code.  As  a  result  we  start  to  have  ambiguous  resolved  directions  of  arrival 
as  observed  in  Figure  5.4.  Here,  the  ratio  XJk=2;  therefore  the  number  of  folds  and 
quantization  levels  are  twice  the  values  found  for  the  ideal  case,  as  illustrated  in  Figure  5.6 
for  modulus  m3=5. 


63 


In  all  cases,  the  number  of  folds  present  in  the  phase  response  waveform  will  be 
proportional  to  the  ratio 

In  summary,  for  the  results  obtained  in  the  third  case,  we  won’t  be  able  to  find  the 
direction  of  arrival  over  the  complete  field  of  view  without  some  ambiguous  results. 
However,  when  the  frequency  of  the  incoming  signal  is  smaller  than  f0,  we  do  not  have 
any  ambiguities,  but  the  overall  error  of  the  resolved  direction  of  arrival  increases. 


-80  -60  -40  -20  O  20  40  60  80 
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Figure  5.6:  Phase  Response  Waveform  (Mod  5)  for  f  =8.5  GHz  (Dashed)  and 

f =17  GHz. 


This  suggests  that,  for  example,  if  one  of  the  requisites  of  our  antenna  is  to  be  able  to  find 
the  direction  of  arrival  of  any  incoming  signal  with  frequencies  over  some  range,  say  2  to 
18  GHz,  the  higher  frequency  in  that  range  should  be  chosen  to  build  our  antenna,  i.e.,  18 
GHz  for  this  case.  By  using  this,  the  direction  of  arrival  information  will  always  be 
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obtained  without  ambiguities  over  the  entire  field  of  view.  However,  the  overall 
quantization  error  will  be  different  for  each  case  depending  on  the  number  of  quantization 
levels  that  will  be  used  and,  consequently,  on  the  frequency  value  of  the  incoming  signal. 

At  this  point,  we  determine  how  to  correct  the  additional  error  when  the  frequency 
of  the  incoming  signal  is  different  than  f0  (frequency  used  to  build  our  antenna).  To  correct 
the  error  however,  the  frequency  of  the  incoming  signal  must  be  known.  With  the 
measured  frequency,  the  number  of  folds  present  in  the  phase  response  waveform  is 
determined.  The  new  number  of  folds  present  in  the  phase  response  over  the  entire  field  of 
view  is  given  by 


A  \ 


ni.fa  = 


n. 


i,RNS* 


(5-3) 


Consider  modulus  3  and  fr=4.25  GHz.  For  this  particular  case,  the  new  number  of  folds 
(ni,fn)  ^  equal  to  10.  From  Equation  (4-3)  and  knowing  that  the  number  of  quantization 
levels  is  equal  to  the  dynamic  range,  the  actual  number  of  quantization  levels  (NQL)  for  a 
specific  frequency  is  given  by 


nql=  ni,f„  ,  (5-4) 

where  n^  is  the  new  number  of  folds  and  n\  represents  the  PRP  modulus  that  is  being 
considered.  For  our  example,  the  number  of  quantization  levels  is  equal  to  30.  With  this 
number  of  quantization  levels  it  is  possible  to  find  the  transition  points  from  one 
quantization  level  to  another.  The  transition  points  will  occur  at  some  angle  0t  for  which  the 
sin(0t)  corresponds  to  some  value  in  the  range  from  0  to  ±1  given  by 
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(5-5) 


x(n)  =  ± - 

N 

^QL 

for  n  =  0, 1,  2,...,  NQL.  Therefore,  the  transition  angles  are  given  by 

0t  (n)  =  sin-1  (x(n)).  (5-6) 

Having  the  transition  angles,  the  midvalue  of  each  quantization  level  can  easily  be  found,  as 
well  as  the  corresponding  quantization  error.  The  new  midvalues  for  each  quantization  level 
are  then  indexed  to  the  RNS  logic  block  so  that  for  each  3  RNS  codes  will  correspond  the 
correct  resolved  direction  of  arrival.  Since  the  phase  response  waveform  is  always  phased 
at  0°  (see  Figures  4.5,  4.6,  and  4.7),  it  is  possible  to  implement  this  new  algorithm  in  our 
RNS  logic  block  to  obtain  the  correct  resolved  direction  of  arrival.  Figure  5.7  shows  the 
resolved  direction  of  arrival  and  the  quantization  error  for  an  incoming  signal  with  a 
frequency  fr=4.25  GHz,  using  this  new  algorithm  to  introduce  the  corrections.  With  this 
new  algorithm,  the  error  in  the  resolved  direction  of  arrival  when  the  frequency  of  the 
incoming  signal  is  smaller  than  f0,  becomes  very  small  when  compared  with  the  results 
illustrated  in  Figure  5.3.  Figure  5.8  shows  the  quantization  error  for  this  case  using  another 
scale. 
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igure  5.7:  Resolved  Direction  of  Arrival  and  Quantization  Error  for  an  Incoming  Signal 

With  a  Frequency  fr=4.25  GHz. 
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Figure  5.8:  Characteristic  Quantization  Error  Curve  for  an  Incoming  Signal  With  a 

Frequency  f=4.25  GHz. 


D.  RELATIONSHIP  BETWEEN  ANTENNA  SIZE,  MODULI 
SET  AND  FREQUENCY  CHOSEN  TO  BUILD  THE 
ANTENNA 

The  size  of  the  antenna  is  an  important  parameter  when  considering  a  new  system  to 
be  implemented.  The  relationship  between  frequency  and  wavelength  is  given  by  A,=c/f 
where  c  is  the  speed  of  light  in  free  space.  This  relationship  is  important  because  the  length 
of  the  DF  antenna  is  given  by  the  distance  between  the  widest  pair  of  elements  which  is 
proportional  to  the  highest  number  of  wavelengths  (see  Equation  (4-8))  required  to 
correctly  apply  the  RNS  code  to  a  certain  PRP  moduli  set.  On  other  hand,  when  building 
the  antenna,  the  highest  frequency  in  the  range  of  frequencies  that  are  being  observed  must 
be  chosen  in  order  to  be  able  to  obtain  the  direction  of  arrival  information  over  the  entire 
field  of  view  without  ambiguities.  Table  5.1  shows  some  different  PRP  moduli  set 
combinations  and  the  distance  in  wavelengths  required  for  the  widest  pair  of  elements  in  the 
array. 


PRP  Moduli  Set 

(m,,  m„  m,) 

Distance  between  the 
widest  pair  of  elements 

(9,  13) 

6.5  X 

(17,  20) 

10  X 

(47,  50) 

25  X 

( 97,  101) 

50.5  X 

(3,  4,  5) 

10  X 

(5,  7,  9) 

31.5  X 

(7,  9,  11) 

49.5  X 

(11,  13,  15) 

97.5  X 

Table  5.1:  Distance  Between  the  Widest  Pair  of  Elements 
According  to  the  Moduli  Set 
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Separation  Between  the  Widest  Pair  of  Elements  (m) 


For  each  case  considered  in  Table  5.1,  Figure  5.9  shows  the  required  distance  between  the 
widest  pair  of  elements  as  a  function  of  the  frequency  used  to  build  the  antenna. 


Figure  5.9:  Required  Distance  Between  the  Widest  Pair  of  Elements  as  a  Function 
of  the  Frequency  Used  to  Build  the  Antenna. 


The  results  show  that  it  is  advantageous  to  chose  a  frequency  as  high  as  possible  to  build 
the  antenna.  First,  with  a  higher  frequency  it  is  possible  to  obtain  the  direction  of  arrival 
information  without  ambiguities  over  a  larger  bandwidth.  In  addition,  for  the  same  moduli 
set,  the  size  of  the  antenna  becomes  smaller  when  using  higher  frequencies  to  build  it. 
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E.  QUANTIZATION  ERROR  BEHAVIOR  VS  MODULI  SET 


Another  interesting  point  is  to  investigate  how  the  quantization  error  varies 
according  to  the  chosen  moduli  set.  The  overall  quantization  error  of  the  system  is  directly 
related  to  the  chosen  PRP  moduli  set,  since  it  determines  the  dynamic  range,  and  therefore 
the  number  of  quantization  levels  available  to  quantize  the  phase  response  over  the  entire 
field  of  view.  With  a  higher  number  of  quantization  levels,  the  quantization  error  is  smaller. 
Therefore,  a  better  performance  of  our  antenna  is  to  be  expected.  Figures  5.10  through 
5.17  show  the  characteristic  quantization  error  curve  associated  with  each  of  the  moduli  set 
considered  in  Table  5.1,  assuming  that  the  antenna  is  build  using  a  frequency  (f0)  equal  to 
18  GHz. 


Input  Direction  of  Arrival  (  degrees  ) 

Figure  5.10:  Characteristic  Quantization  Error  Curve  for  the  PRP  Moduli  Set  (9, 13). 
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Figure  5.1 1:  Characteristic  Quantization  Error  Curve  for  the  PRP  Moduli  Set  (17, 20). 
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Figure  5.12:  Characteristic  Quantization  Error  Curve  for  the  PRP  Moduli  Set  (47, 50). 
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Figure  5.13:  Characteristic  Quantization  Error  Curve  for  the  PRP  Moduli  Set  (97, 101). 
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Figure  5.15:  Characteristic  Quantization  Error  Curve  for  the  PRP  Moduli  Set  (5, 7,  9). 
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Figure  5.17:  Characteristic  Quantization  Error  Curve  for  the  PRP  Moduli  Set  (1 1, 13, 15). 


Comparing  the  results  obtained  for  each  case  it  is  clear  that  when  using  a  N  moduli 
set  with  N=2,  the  quantization  error  is  smaller  than  when  using  a  moduli  set  with  N=3, 
assuming  that  the  widest  distance  between  elements  is  the  same.  Therefore,  with  a  moduli 
set  with  N=2,  it  is  possible  to  obtain  a  very  high  accuracy  with  relative  small  antenna  sizes. 
However,  the  required  number  of  comparators  for  this  scheme  (N=2)  becomes  very  high 
when  compared  with  the  scheme  using  a  moduli  set  with  N=3.  Table  5.2  contains  a 
compilation  of  all  information  of  each  case  considered  in  Table  5.1.  In  addition,  the  last 
column  shows  the  equivalent  rms  accuracy  of  an  interferometer  (conventional  methods) 
assuming  that  the  direction  of  arrival  is  equal  to  0°  and  the  SNR  is  equal  to  0  dB.  The  rms 
angular  accuracy  of  an  interferometer  in  radians  [Ref.  1]  is  given  by 


_  Aa 
°9  “  tWSNR  ’ 


(5-7) 
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where  Aa  is  the  separation  between  adjacent  nulls  given  by 


Aa  = 


X 

dcosG 


(5-8) 


Moduli  Set 

Largest 

dist. 

— w. — 

Comparators 

Min.  Quant. 
Error  (deg.) 
fr=  f0=18  GHz 

Min.  Quant. 
Error  (deg.) 
fr  =  2  GHz 

(9, 13) 

6.5  X 

20 

±  0.4897 

±  4.425 

4.897 

(17, 20) 

10  X 

35 

±0.1685 

±1.517 

1.827 

(47, 50) 

25  X 

95 

±  0.0244 

±  0.2194 

0.729 

(97, 101) 

50.5  X 

196 

±  0.0058 

±  0.0526 

0.365 

(3,4,5) 

10  X 

9 

±  0.9551 

±  8.729 

1.827 

Kftftn 

31.5  X 

18 

±0.1819 

±  1.638 

0.579 

(7,9,11) 

49.5  X 

24 

±  0.0827 

±  0.744 

0.368 

(11,13,15) 

97.5  X 

36 

±  0.0267 

±  0.2404 

0.187 

Table  5.2:  Comparison  Between  Different  RNS  Architecture  Schemes  and  the 

Conventional  Methods. 


As  shown  in  Table  5.2,  with  this  new  architecture  it  is  possible  to  obtain  the  direction  of 
arrival  information  with  a  substantial  improvement  in  accuracy  when  compared  with  the 
conventional  interferometer  methods. 


75 


76 


VI.  SUMMARY  AND  CONCLUDING  REMARKS 


In  this  thesis  two  novel  architectures  for  high-resolution  wideband  direction  finding 
based  in  the  SNS  and  RNS  were  developed  and  investigated.  The  first  architecture 
incorporates  a  SNS  encoding  of  the  amplitude  of  the  electrical  field  response  (symmetrical 
folding  waveform)  of  the  phase  response  (sawtooth  folding  waveform)  array.  The  second 
architecture  incorporates  a  RNS  encoding  of  the  interferometer  phase.  In  addition,  an  array 
antenna  was  constructed  in  a  ground  plane  using  rectangular  waveguide  elements  with  a 
center  frequency  f0=8.5  GHz  and  the  SNS  architecture  was  tested. 

The  simulated  results  show  that  the  SNS  architecture  is  sensitive  to  the  drop  off  in 
the  antenna  gain  due  to  the  element  factor  effect  associated  with  the  kind  of  individual 
elements  used  to  build  the  antenna.  In  this  case  it  is  not  possible  to  obtain  the  direction  of 
arrival  information  over  the  entire  field  of  view  (from  0°  to  ±90°).  However,  we  still  have  a 
small  interval  (from  0°  to  ±20°)  where  it  is  possible  to  have  that  information.  The 
experimental  results  indicate  a  considerable  agreement  for  this  case. 

For  the  RNS  architecture  the  simulated  results  indicate  that  it  provides  a  more 
robust  solution  to  high-resolution  DF  having  an  excellent  performance.  In  fact,  since  the 
phase  response  has  an  uniform”  distribution  it  is  possible  to  obtain  the  direction  of  arrival 
information  over  the  entire  field  of  view.  To  avoid  ambiguities  in  the  resolved  direction  of 
arrival  (over  the  entire  field  of  view)  the  spacing  between  elements  must  be  established  for 
the  highest  frequency  to  be  received.  On  other  hand,  the  frequency  of  the  incoming  signal 
has  to  be  known  in  order  to  make  the  necessary  corrections  in  the  RNS  logic  block  using  a 
fast  correction  algorithm  to  obtain  the  correct  direction  of  arrival.  The  performance  of  the 
algorithm  was  quantified  and  shown  to  have  excellent  results  over  a  large  bandwidth. 
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The  quantization  error  associated  with  both  architectures  depends  on  the  moduli  set 
that  is  chosen  to  build  the  antenna,  which  also  determines  the  size  of  the  antenna  and  the 
required  number  of  comparators.  Table  5.2,  in  Chapter  V,  illustrates  some  specifications 
for  the  RNS  architecture  using  different  moduli  set  combinations.  It  also  shows  the 
equivalent  rms  accuracy  for  the  conventional  methods  according  to  the  distance  between  the 
widest  pair  of  elements.  Comparing  those  values  with  the  values  found  for  the  RNS 
architecture  we  notice  a  great  improvement  in  accuracy  with  the  new  architecture. 

The  advantages  of  the  RNS  architecture  over  the  conventional  methods  can  be 
summarized  as  follows: 

•  Higher  resolution,  higher  bandwidth  capability, 

•  Small  number  of  antenna  elements, 

•  Small  array  size  to  achieve  higher  accuracy, 

•  Only  few  comparators  are  required,  depending  on  the  chosen  moduli  set, 

•  Full  coverage  from  0°  to  ±90°  without  ambiguities. 

It  was  not  possible  to  obtain  acceptable  experimental  results  for  the  RNS  antenna 
architecture.  The  main  problem  was  the  nonavailability  of  the  phase  detectors.  However, 
future  investigations  are  currently  underway  and  will  correct  this  and  another  problems. 
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APPENDIX:  MATLAB™  CODE 


%Program  1 

%This  program  'Ampanal.m'  calculates  the  resolved  Direction  of  Arrival  using  the  SNS 
%code  to  process  the  amplitude  response  of  the  array  for  the  moduli  set  (3, 4, 5). 

clear  all 
whitebg 
global  ml 
global  m2 
global  m3 
global  fo 
global  fr 
global  ds 
global  rad 
global  M 
global  theta 

obj=input('Do  you  want  the  Array  Factor  (1)  or  the  Electrical  field  (2)  analysis  ?') 
ds=input('What  is  the  step  size?  ') 

fo=input('What  was  the  frequency  used  to  built  the  antenna  (GHz)  ?') 
fr=inputCWhat  is  the  frequency  of  the  received  signal  (GHz)  ?') 


ml  =3; 
m2  =4; 
m3=5; 
fo=fo*le9; 
fr=fr*le9; 


rad=pi/180; 

theta=-90:ds:90; 

c=3e8; 

L=c/fo; 

Ll=c/fr; 


M=ml*m2*m3; 


^********Hc**********£)jstances  Between  Each  Pair  of  Elements****************** 

dl=0.5*M/ml*L; 

d2=0.5*M/m2*L; 

d3=0.5*M/m3*L; 


%*****Generates  the  Array  Factor  or  the  Electrical  field  for  Each  Pair  of  Elements****** 

phil=(2*pi*dl/L)*sin(theta*rad); 

phi2=(2*pi*d2/L)  *  sin(theta*rad); 

phi3=(2*pi*d3/L)*sin(theta*rad); 


if  obj==l 

modlA=cos(0.5*(phil)+pi); 
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mod2A=cos(0.5*(phi2)+pi/2); 

mod3A=cos(0.5*(phi3)+pi); 

else 

modlA=cos(0.5*(phil)+pi).*cos(theta*rad); 

mod2A=cos(0.5*(phi2)+pi/2).*cos(theta*rad); 

mod3A=cos(0.5*(phi3)+pi).*cos(theta*rad); 

end 

figure(l) 
subplot(311) 
plot(theta,mod  1  A) 
axis([-90,90,-l,l]); 

xlabel('Input  Direction  of  Arrival  ( degrees )') 
grid 

subplot(312) 
plot(theta,mod2A) 
axis([-90,90,- 1,1]); 

xlabel('Input  Direction  of  Arrival  ( degrees )') 

ylabel(' Array  Factor  Amplitude ') 

grid 

subplot(313) 

plot(theta,mod3A) 

axis([-90,90,-l,l]); 

xlabel(’Input  Direction  of  Arrival  ( degrees )') 
grid 

^)*********s(!****(];ajcu2a^es  Threshold  Values  for  Each  Modulus**************** 
thresa 

%*********Qenerates  the  Comparators  RNS  Output  code  for  Modulus  rnl=3********** 
for  i=l:180/ds+l 
if  modi  A(i)<=Tl(l) 
modl(i)=0; 

elseif  modlA(i)>Tl(2) 

modl(i)=2; 

else 

modl(i)=l; 

end 

end 

end 

%******** *Generates  the  Comparators  RNS  Output  code  for  Modulus  ni2=4********** 
for  i=l:180/ds+l 
if  mod2A(i)<=T2(l) 
mod2(i)=0; 

elseif  mod2A(i)<=T2(2) 
mod2(i)=l; 

elseif  mod2A(i)<=T2(3) 
mod2(i)=2; 
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else 

mod2(i)=3; 

end 

end 

%*********Generates  the  Comparators  RNS  Output  code  for  Modulus  rn3=5********** 
for  i=l:180/ds+l 
if  mod3A(i)<=T3(l) 
mod3(i)=0; 

elseif  mod3A(i)<=T3(2) 
mod3(i)=l; 

elseif  mod3A(i)<=T3(3) 
mod3(i)=2; 

elseif  mod3A(i)<=T3(4) 

mod3(i)=3; 

else 

mod3(i)=4; 

end 

end 

%****  *  ^Calculates  the  Resolved  Direction  of  Arrival  for  the  moduli  Set  (3, 4, 5)******** 
SNSLogicBlock 

%*Plots  the  Resolved  Direction  of  Arrival  as  a  Function  of  the  Input  Direction  of  Arrival** 

figure(2) 

plot(theta,th,'g') 

xlabel(Input  Direction  of  Arrival  ( degrees )') 
ylabel('Resolved  Direction  of  Arrival  ( degrees )') 
grid 

axis([-90  90  -90  90]) 
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%  Program  2 

%  This  program  ‘thresa.m’  calculates  the  Threshold  Values  for  the  moduli  set  (3,  4,  5) 
%when  analyzing  the  Amplitude  Response. 

global  ml 
global  m2 
global  m3 
global  fo 
global  ds 
global  M 
global  L 


%*******CalCUlateS  the  transition  angles  from  one  to  another  quantization  level  ******** 

x=-l:2/M:l; 

y=asin(x); 


%**************Cajcula|.es  Threshold  Values  for  modulus  1x11=3*************** 

for  l=l:ml-l 

angl(l)=y(l+l); 

phil  (l)=(2*pi*dl/L)*sin(angl(l)); 

T 1  (l)=cos(0.5*(phil  (l))+pi); 
end 

c^*******^:^^****^^^  hxe  Threshold  Values  for  Modulus  m2=4*************** 

forr=l:m2-l 

ang2(r)=y(r+l); 

phi2(r)=(2*pi*d2/L)*sin(ang2(r)); 

T2(r)=cos(0.5*(phi2(r))+pi/2); 

end 

%:i:*************caicuiateS  the  Threshold  Values  for  Modulus  m3=5*************** 

for  h=l:m3-l 

ang3(h)=y(h+l); 

phi3(h)=(2*pi*d3/L)*sin(ang3(h)); 

T3(h)=cos(0.5*(phi3(h))+pi); 

end 
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%Program  3 

%This  program  ‘SNSLogicBlock.m’  calculates  the  resolved  direction  of  arrival  according 
to  the  3-SNS  code. 


global  ds 
global  ml 
global  m2 
global  m3 
global  fo 
global  fr 
global  M 
global  theta 

%****Cainbines  the  SNS  Code  coming  from  each  channel  ml=3,  m2=4,  and  m3=5***** 
res=[modl'  mod2'  mod3']; 

%***Calculates  the  Resolved  Direction  of  Arrival  according  to  the  3-SNS  Code  value**** 
for  u=l:180/ds+l 
if  res(u,:)==([0,0,0]) 
th(u)=-82.5824; 
elseif  res(u,:)= ([1,1,1]) 
th(u)=-72.0627; 
elseif  res(u,:)==([2,2,2]) 
th(u)=-66.5593; 
elseif  res(u,:)=([2,3,3]) 
th(u)=-62.1158; 
elseif  res(u, :)==([  1,3 ,4]) 
th(u)=-58.2581; 
elseif  res(u,:)==([0,2,4]) 
th(u)=-54.7864; 
elseif  res(u,:)==([0,l,3]) 
th(u)=-5 1.5928; 
elseif  res(u,:)==([l,0,2]) 
th(u)=-48.6110; 
elseif  res(u,:)==([2,0,l]) 
th(u)=-45.7968; 
elseif  res(u,:)==([2,l,0]) 
th(u)=-43.1187; 
elseif  res(u,:)==([l,2,0]) 
th(u)=-40.5534; 
elseif  res(u,:)=([0,3,l]) 
th(u)=-38.0832; 
elseif  res(u,:)==([0,3,2]) 
th(u)=-35.6940; 
elseif  res(u,:)=([l,2,3]) 
th(u)=-33.3745; 
elseif  res(u,:)==([2,l,4]) 
th(u)=-31.1155; 
elseif  res(u,:)==([2,0,4]) 
th(u)=-28.9091; 
elseif  res(u,  :)==([!,  0,3]) 
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th(u)=-26.7487; 

elseif  res(u,:)=([0,l,2]) 

th(u)=-24.6287; 

elseif  res(u,:)==([0,2,l]) 

th(u)=-22.5442; 

elseif  res(u,:)==([l,3,0]) 

th(u)=- 20.4907; 

elseif  res(u,:)==([2,3,0]) 

th(u)=- 18.4644; 

elseif  res(u,:)==([2,2,l]) 

th(u)=-16.4618; 

elseif  res(u,:)=([l,l,2]) 

th(u)=- 14.4797; 

elseif  res(u,;)==([0,0,3]) 

th(u)=-12.5152; 

elseif  res(u,:)==([0,0,4]) 

th(u)=- 10.5655; 

elseif  res(u,;)=([l,l,4]) 

th(u)=-8.6382; 

elseif  res(u,:)==([2,2,3]) 

th(u)=-6.7007; 

elseif  res(u,:)==([2,3,2]) 

th(u)=-4.7809; 

elseif  res(u,:)==([l,3,l]) 

th(u)=-2.8664; 

elseif  res(u,:)==([0,2,0]) 

th(u)=-0.9551; 

elseif  res(u,:)==([0,l,0]) 

th(u)=0.9551; 

elseif  res(u,;)==([l,0,l]) 

th(u)=2.8664; 

elseif  res(u,:)==([2,0,2]) 

th(u)=4.7809; 

elseif  res(u,:)==([2,l,3]) 

th(u)=6.7007; 

elseif  res(u,:)==([l,2,4]) 

th(u)=8.6382; 

elseif  res(u,:)==([0,3,4]) 

th(u)=10.5655; 

elseif  res(u,;)==([0,3,3]) 

th(u)=12.5152; 

elseif  res(u,:)==([l,2,2]) 

th(u)=14.4797; 

elseif  res(u,;)==([2,l,l]) 

th(u)=16.4618; 

elseif  res(u,:)==([2,0,0]) 


th(u)=l  8.4644; 


elseif  res(u,:)==([l,0,0]) 
th(u)=20.4907; 


elseif  res(u,:)==([0,l,l]) 
th(u)=22.5442; 


I 
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elseif  res(u,:)==([0,2,2]) 

th(u)=24.6287; 

elseif  res(u,:)=([l,3,3]) 

th(u)=26.7487; 

elseif  res(u,:)==([2,3,4]) 

th(u)=28.9091; 

elseif  res(u,:)==([2,2,4]) 

th(u)=31.1155; 

elseif  res(u,:)==([l,l,3]) 

th(u)=33.3745; 

elseif  res(u,:)=([0,0,2]) 

th(u)=35.6940; 

elseif  res(u,:)==([0,0,l]) 

th(u)=38.0832; 

elseif  res(u,:)=([l, 1,0]) 

th(u)=40.5534; 

elseif  res(u,:)==([2,2,0]) 

th(u)=43.1187; 

elseif  res(u,:)==([2,3,l]) 

th(u)=45.7968; 

elseif  res (u, :)==([  1,3, 2]) 

th(u)=48.6110; 

elseif  res(u,:)==([0,2,3]) 

th(u)=51.5928; 

elseif  res(u,:)=([0,l,4]) 

th(u)=54.8;%54; 

elseif  res(u,:)==([l,0,4]) 

th(u)=58.2581; 

elseif  res(u,:)==([2,0,3]) 

th(u)=62.1158; 

elseif  res(u,:)==([2,l,2]) 

th(u)=66.5593; 

elseif  res(u,:)==([  1,2,1]) 

th(u)=72.0627; 

elseif  res(u,:)=([0,3,0]) 

th(u)=82.5824; 

end 

end 
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%  Program  4 

%This  program  'Phanalnew'  calculates  the  Resolved  Direction  of  Arrival  using  the  RNS 
%code  to  process  the  phase  response  of  the  array  for  the  moduli  set  (3, 4,  5). 

clear  all 
whitebg 
global  ml 
global  m2 
global  m3 
global  fo 
global  fr 
global  ds 
global  rad 
global  M 
global  theta 

ds=input('What  is  the  step  size?  ') 

fo=input('What  was  the  frequency  used  to  built  the  antenna  (GHz)  ?') 
fr=inputCWhat  is  the  frequency  of  the  received  signal  (GHz)  ?') 


ml=3; 

m2=4; 

m3=5; 

fo=fo*le9; 

fr=ff*le9; 

rad=pi/180; 

theta=-90:ds:90; 

c=3e8; 

L=c/fo; 

Ll=c/fr; 


M=ml*m2*m3; 


%*****=c*******^*****Distarices  getween  Each  Pair  of  Elements****************** 

dl=0.5*M/ml*L; 

d2=0.5*M/m2*L; 

d3=0.5*M/m3*L; 


«^**************Qenerates  An-ay  Factor  for  Each  Pair  of  Elements************** 

phil=(2*pi*dl/Ll)*sin(theta*rad); 

modlp=angle(exp(j*((phil+pi)))); 

phil=(2*pi*d2/Ll)*sin(theta*rad); 
mod2p=angle(exp(j  *((phi  1  )+pi))); 

phil=(2*pi*d3/Ll)*sin(theta*rad); 
mod3p=angle(exp(j  *((phi  1  )+pi))); 
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figure(l) 

subplot(311) 

plot(theta,modlp) 

axis([-90,90,-4,4]); 

xlabel('Input  Direction  of  Arrival  ( degrees )') 
grid 

subplot(312) 

plot(theta,mod2p) 

axis([-90,90,-4,4]); 

xlabel('Input  Direction  of  Arrival  ( degrees )') 

ylabel('Phase  Response  ( pi ) ') 

grid 

subplot(313) 

plot(thetaimod3p) 

axis([-90,90,-4,4]); 

xlabel('Input  Direction  of  Arrival  ( degrees )') 
grid 

%**************Qajcujates  the  Threshold  Values  for  Each  Modulus**************** 
thres 

%*********Generates  the  Comparators  RNS  Output  code  for  Modulus  rnl=3********** 
for  i=l:180/ds+l 
if  modlp(i)<=Tl(l) 
modl(i)=0; 
elseif  modlp(i)>Tl(2) 
modl(i)=2; 
else 

modl(i)=l; 

end 

end 

end 

%*********Generates  the  Comparators  RNS  Output  code  for  Modulus  rn2=4********** 
for  i=l:180/ds+l 
if  mod2p(i)<=T2(l) 
mod2(i)=0; 

elseif  mod2p(i)<=T2(2) 
mod2(i)=l; 

elseif  mod2p(i)<=T2(3) 

mod2(i)=2; 

else 

mod2(i)=3; 

end 

end 

%*********Generates  the  Comparators  RNS  Output  code  for  Modulus  ni3=5********** 
for  i=l:180/ds+l 
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if  mod3p(i)<=T3(l) 
mod3(i)=0; 

elseif  mod3p(i)<=T3(2) 
mod3(i)=l; 

elseif  mod3p(i)<=T3(3) 
mod3(i)=2; 

elseif  mod3p(i)<=T3(4) 

mod3(i)=3; 

else 

mod3(i)=4; 

end 

end 

%******CalCUlate  the  Resolved  Direction  of  Arrival  and  the  Quantization  Error********* 

^*****!(:^******!i:*H:********£or  ^  get  ^  ^  ^************************* 

RNSlogicBlock 

%*********P20ts  ^  Resoive(j  Direction  of  Arrival  and  the  Quantization  Error********* 

%*********=M=******as  a  function  of  the  Input  Direction  of  Arrival***************** 

figure(2) 

plot(theta,th,'g') 

hold 

xlabel(’Input  Direction  of  Arrival  ( degrees )') 

ylabel('Resolved  DOA  and  Quantization  Error  ( degrees )') 

axis([-90  90  -90  90]) 

plot(y([2:size(y,2)]),k,'r') 

gnd 
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%  Program  5 

%This  program  ‘thres.m’  calculates  the  Threshold  values  for  the  moduli  set  (3,  4,  5) 
%when  analyzing  the  phase  response. 

global  ml 
global  m2 
global  m3 
global  fo 
global  fir 
global  ds 
global  M 
global  L 


%*******CalCUlate  the  transition  angles  from  one  to  another  quantization  level********* 

x=-l:2/M:l; 

y=asin(x); 


%%*:):***********Qaicuiate  the  Threshold  Values  for  modulus  ml=3**************** 

for  l=l:ml-l 

angl(l)=y(l+l); 

phil(l)=(2*pi*dl/L)*sin(angl(l)); 

T1  (l)=angle(exp(j  *((phi  1  (l))+pi))); 
end 

^**************(231^!!^  the  Threshold  Values  for  modulus  m2=4**************** 

for  r=l:m2-l 

ang2(r)=y(r+l); 

phi2(r)=(2*pi*d2/L)*sin(ang2(r)); 

T2(r)=angle(exp(j*((phi2(r))+pi))); 

end 


^******>1!*******^^^  t^e  Threshold  Values  for  modulus  m3=5**************** 

forh=l:m3-l 

ang3(h)=y(h+l); 

phi3(h)=(2*pi*d3/L)*sin(ang3(h)); 

T3(h)=angle(exp(j*((phi3(h))+pi))); 

end 
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%  Program  6 

%This  program  ‘RNSlogicBlock.m’  calculates  the  Resolved  Direction  of  Arrival  and 
%Quantization  Error  according  to  the  3-RNS  code.  It  also  uses  an  algorithm  to  calculate  the 
%correct  quantization  level  midvalue  when  the  frequency  of  the  received  signal  (fr)  is 
%different  form  the  center  frequency  (fo)  used  in  the  antenna  design. 


global  rad 
global  ml 
global  m2 
global  m3 
global  fo 
global  fr 
global  M 
global  theta 


%*****CalCulate  the  total  number  of  quantization  levels  over  the  entire  field  of  view***** 

%***********an(j  (jje  transition  angles  from  one  to  another  quantization  level********** 

nf=fr/fo*M/ml; 

nql=nf*ml; 

p=0:2/nql:l; 

if  max(p)~=l 

P=[p  l]; 

end 

p=[  -rot90(p(2:size(p,2)))'  p]; 
trang=  l/rad*asin(p); 

%** Calculate  the  quantization  level  width  and  the  midvalues  for  each  quantization  level*** 
for  n=l:size(p,2)-l 

midval(n)=trang(n)+(trang(n+l)-trang(n))/2; 

qlw(n)=midval(n)-trang(n); 

end 

if  rem(size(midval,2),2)==0 

a=midval; 

b=size(a,2)/2; 

end 

if  rem(size(midval,2),2)==l 

a=midval; 

b=(size(a,2)+l)/2; 

end 

if  size(a,2)~=60 
a=[zeros(l,b)  a  zeros(l,b)]; 
end 

b=size(a,2)/2; 

for  n=2:2*size(trang,2)-l 

if  rem(n,2)==0 

y(n)=trang((n)/2); 

else 

y  (n)=trang((n+ 1  )/2); 
end 
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end 


i^*********************^^^^  Quantization  Error*********************** 

for  n=l:size(trang,2)-l 

c(n)=(trang(n+l)-trang(n))/2; 

end 

for  r=l:2*size(trang,2)-2 
if  rem(r,2)==0 
k(r)=c(r/2); 

k(r)=-c((r+l)/2); 

end 

end 


%****Combines  the  RNS  Code  coming  from  each  channel  ml=3,  m2=4,  and  m3=5***** 

res=[mod3'  mod4'  mod5']; 

for  u=l:180/ds+l 
if  res(u,:)==([0,0,0]) 
th(u)=a(b+l); 
elseif  res(u,:)==([l, 1,1]) 
th(u)=a(b+2); 
elseif  res(u,:)=([2,2,2]) 
th(u)=a(b+3); 
elseif  res(u,:)— ([0,3,3]) 
th(u)=a(b+4); 
elseif  res(u,:)=([l,0,4]) 
th(u)=a(b+5); 
elseif  res(u,:)==([2, 1,0]) 
th(u)=a(b+6); 
elseif  res(u,:)==([0,2,l]) 
th(u)=a(b+7); 
elseif  res(u,:)==([l,3, 2]) 
th(u)=a(b+8); 
elseif  res(u,:)=([2,0,3]) 
th(u)=a(b+9); 
elseif  res(u,:)=([0, 1,4]) 
th(u)=a(b+10); 
elseif  res(u,:)==([l,2,0]) 
th(u)=a(b+ll); 
elseif  res(u,:)==([2,3,l]) 
th(u)=a(b+12); 
elseif  res(u,:)==([0,0,2]) 
th(u)=a(b+13); 

elseif  res(u,:)==([l,l, 3]) 

th(u)=a(b+14); 

elseif  res(u,:)==([2,2,4]) 

th(u)=a(b+15); 

elseif  res(u,:)==([0,3,0]) 

th(u)=a(b+16); 

elseif  res(u,:)==([l,0,l]) 
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th(u)=a(b+17); 

elseif  res(u,:)==([2,l,2]) 

th(u)=a(b+18); 

elseif  res(u,:)==([0,2,3]) 

th(u)=a(b+19); 

elseif  res(u,:)==([l,3,4]) 

th(u)=a(b+20); 

elseif  res(u,:)==([2,0,0]) 

th(u)=a(b+21); 

elseif  res(u,:)==([0,l,l]) 

th(u)=a(b+22); 

elseif  res(u,:)==([l,2,2]) 

th(u)=a(b+23); 

elseif  res(u,:)==([2,3,3]) 

th(u)=a(b+24); 

elseif  res(u,:)==([0,0,4]) 

th(u)=a(b+25); 

elseif  res(u, :)==([  1,1,0]) 

th(u)=a(b+26); 

elseif  res(u,:)==([2,2,l]) 

th(u)=a(b+27); 

elseif  res(u,:)==([0,3,2]) 

th(u)=a(b+28); 

elseif  res(u,:)==([l,0,3]) 

th(u)=a(b+29); 

elseif  res(u,:)==([2,l,4]) 

th(u)=a(b+30); 

elseif  res(u,:)==([0,2,0]) 

th(u)=a(b-29); 

elseif  res(u,:)==([l,3,l]) 

th(u)=a(b-28); 

elseif  res(u,:)==([2,0,2]) 

th(u)=a(b-27); 

elseif  res(u,:)==([0,l,3]) 

th(u)=a(b-26); 

elseif  res(u, :)==([!, 2, 4]) 


th(u)=a(b-25); 

elseif  res(u,:)==([2,3,0]) 

th(u)=a(b-24); 


elseif  res(u,:)==([0,0,l]) 
th(u)=a(b-23); 


elseif  res(u,:)==([l,l,2]) 
th(u)=a(b-22); 


elseif  res(u,:)==([2,2,3]) 
th(u)=a(b-21); 


elseif  res(u,:)==([0,3,4]) 
th(u)=a(b-20); 
elseif  res(u,:)==([l,0,0]) 
th(u)=a(b-19); 


elseif  res(u,:)==([2,l,l]) 
th(u)=a(b-18); 


elseif  res(u,:)==([0,2,2]) 
th(u)=a(b-17); 
elseif  res(u, :)==([  1,3, 3]) 
th(u)=a(b-16); 
elseif  res(u,:)==([2,0,4]) 
th(u)=a(b-15); 
elseif  res(u,:)==([0,l,0]) 
th(u)=a(b-14); 
elseif  res(u,:)==([l,2,l]) 
th(u)=a(b-13); 
elseif  res(u,:)==([2,3,2]) 
th(u)=a(b-12); 
elseif  res(u,:)==([0,0,3]) 
th(u)=a(b-ll); 
elseif  res(u,:)==([l,l,4]) 
th(u)=a(b-10); 
elseif  res(u,:)==([2,2,0]) 
th(u)=a(b-9); 
elseif  res(u,:)==([0,3,l]) 
th(u)=a(b-8); 
elseif  res(u,:)==([l,0,2]) 
th(u)=a(b-7); 
elseif  res(u,:)==([2,l,3]) 
th(u)=a(b-6); 
elseif  res(u,:)==([0,2,4]) 
th(u)=a(b-5); 
elseif  res(u,:)==([l,3,0]) 
th(u)=a(b-4); 
elseif  res(u,:)==([2,0,l]) 
th(u)=a(b-3); 
elseif  res(u,:)=([0,l,2]) 
th(u)=a(b-2); 
elseif  res(u,:)=([l,2,3]) 
th(u)=a(b-l); 
elseif  res(u,:)==([2,3,4]) 
th(u)=a(b); 
end 
end 
end 
end 


%Program  7 

%This  program  'Quanterror'  calculates  the  Quantization  Error  for  any  PRP  N  moduli  set, 
%  where  N  is  equal  to  2  or  3,  and  any  frequency. 

clear  all 
whitebg 

ds=input('What  is  the  step  size?  ') 

N=input('What  is  the  moduli  set  size  (2  or  3)  ?') 
if  N==2 

ml=input('What  is  the  value  of  modulus  ml  ?') 
m2=input('What  is  the  value  of  modulus  m2  ?') 
end 

if  N==3 

ml=input('What  is  the  value  of  modulus  ml  ?') 
m2=input('What  is  the  value  of  modulus  m2  ?’) 
m3=input('What  is  the  value  of  modulus  m3  ?') 
end 

fo=input('What  was  the  frequency  used  to  built  the  antenna  (GHz)  ?') 
fr=inputeWhat  is  the  frequency  of  the  received  signal  (GHz)  ?') 

%**********NUmber  of  Folds  Over  the  Entire  field  of  View  for  Modulus  ml********** 
nf=fr/fo*M/ml; 

%********Totai  Number  of  Quantization  Levels  Over  the  Entire  Field  of  View********* 
nql=nf*ml; 


%************>i!Generates  the  Quantization  Level  Transistion  Angles**************** 

p=0:2/nql:l; 

if  max(p)~=l 

p=[p  i]; 

end 

p=[  -rot90(p(2:size(p,2)))'  p]; 
trang=l/rad*asin(p); 

^lx:****:f:***!|!****Gaicu2ates  the  Midvalue  of  Each  Quantization  Level**************** 
for  n=l:size(p,2)-l 

midval(n)=trang(n)+(trang(n+l)-trang(n))/2; 

qlw(n)=midval(n)-trang(n) ; 

end 

%********************Gahajiates  the  Quantization  Error************************ 

for  n=2:2*size(trang,2)-l 

if  rem(n,2)==0 

y(n)=trang((n)/2); 

else 

y  (n)=trang((n+ 1  )/2); 

end 

end 
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for  n=l:size(trang,2)-l 
c(n)=(trang(n+ 1  )-trang(n))/2; 
end 

for  r=l:2*size(trang,2)-2 
if  rem(r,2)==0 
k(r)=c(r/2); 

k(r)=-c((r+l)/2); 

end 

end 

%**********plots  the  Resoived  Direction  of  arrival  and  Quantization  Error************ 
%**=M=************as  a  Functjon  0f  the  Input  Direction  of  Arrival****************** 
figure(l) 

pl°t(y([2:size(y,2)]),k,'r') 

xlabel('Input  Direction  of  Arrival  ( degrees )') 

ylabel('  Quantization  Error  ( degrees )') 

axis([-90  90  -round(max(k)+0.5)  round(max(k)+0.5)]) 

grid 
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